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Use  of  Non-SI  Units  in  a Non-NIST  Publication 


It  is  the  policy  of  the  National  Institute  of  Standards  and  Technology  to  use  the  International 
System  of  Units  (metric  units)  in  all  of  its  publications.  However,  in  North  America  in  the 
HVAC&R  industry,  certain  non-SI  units  are  so  widely  used  instead  of  SI  units  that  it  is  more 
practical  and  less  confusing  to  use  measurement  values  for  customary  units  only  in  figures  and 
tables  describing  system  performance. 


EXECUTIVE  SUMMARY 


The  main  goal  of  this  study  was  to  investigate  the  benefits  possible  for  finned  tube  refrigerant 
evaporators  when  refrigerant  distribution  was  precisely  controlled  to  produce  a desired  equal 
superheat  in  each  circuit.  This  goal  was  accomplished  by  examining  three  different  finned  tube 
evaporators;  a wavy  fin,  wavy-lanced  fin,  and  a wavy-lanced  fin  evaporator  with  tube  sheets 
separated.  The  effects  of  non-uniform  airflow  on  capacity  were  also  examined  while  superheat 
was  controlled  in  each  evaporator  circuit.  In  parallel  with  the  experimental  effort,  a modeling 
program  was  implemented  and  validated  with  the  experimental  results  and  then  used  to 
determine  the  savings  in  evaporator  core  volume  possible  if  refrigerant  distribution  was 
controlled  by  a smart  distributor.  In  extreme  cases,  the  savings  in  core  volume  could  be  as  much 
as  40  %. 

Within  the  experimental  part  of  this  study,  all  three  evaporators  could  avoid  significant 
performance  degradation  using  the  ability  to  control  superheat  within  each  of  the  three  finned 
tube  circuits.  As  an  example,  with  cross-counter  flow  configuration,  uniform  airflow,  and  exit 
manifold  superheat  fixed  at  5.6  °C  (10.0  °F),  the  wavy  fin  and  wavy-lanced  fin  evaporator’s 
capacity  dropped  by  as  much  as  41  % and  32  %,  respectively,  as  the  superheat  was  allowed  to 
vary  between  the  circuits.  Control  of  superheat  was  shown  to  be  even  more  important  during 
cross-parallel  refrigerant  flow  due  to  the  rapid  pinching  of  the  refrigerant  and  air  temperatures. 
For  the  wavy  and  lanced  finned  evaporators  in  cross-parallel  flow,  capacity  dropped  by  85  % and 
78  % as  superheat  changed  from  5.6  °C  (10.0  °F)  to  16.7  °C  (30.0  °F).  As  the  coil  faces  were 
blocked  to  produce  a non-uniform  airflow,  pressure  drop  through  the  coils  increased 
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substantially  and  control  of  superheat  was  shown  to  restore  performance.  The  non-uniform 
airflow  tests  showed  that  when  airflow  rate  was  held  constant,  the  losses  in  capacity  due  to  low 
airflow  over  a portion  of  the  coil  could  be  recovered  to  within  2%  of  the  original  uniform  airflow 
capacity  by  controlling  superheat.  The  more  non-uniform  the  airflow  over  the  coil,  the  more 
capacity  was  improved  by  controlling  superheat. 

A combination  of  results  obtained  from  laboratory  testing  and  simulations  indicate  the  influence 
of  tube-to-tube  heat  transfer  on  capacity  degradation.  The  impact  of  tube-to-tube  heat  transfer 
was  negligible  in  tests  with  a uniform  5.6  °C  (10  °F)  superheat,  but  it  was  significant  in  tests 
involving  16.7  °C  (30  °F)  superheat.  Between  the  two  possible  conduction  mechanisms  of  heat 
transfer  that  may  occur,  longitudinal  fin  conduction  is  responsible  for  degraded  performance 
rather  than  longitudinal  tube  conduction,  which  has  insignificant  impact.  The  upgraded  version 
of  the  EVAP5  evaporator  model,  which  accounts  for  tube-to-tube  heat  transfer  based  on  tube 
temperatures,  was  able  to  predict  key  return  bend  temperatures  which  indicated  the  occurrence  of 
tube-to-tube  heat  transfer.  However,  the  study  also  confirmed  that  longitudinal  heat  conduction 
is  affected  by  the  fin  design,  air-side  heat  transfer  coefficient,  and  moisture  removal  process. 
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NOMENCLATURE 


9 9 

Af  = finned  surface  area,  m (ft ) 

A0  = outside  tube  and  fin  area,  m2  (ft2) 

APm  = pipe  mean  surface  area,  m (ft  ) 

Apo  = pipe  outside  surface  area,  m (ft  ) 

Aw  = cross  sectional  area  of  the  tube  available  for  axial  heat  conduction,  m2  (ft2) 

Cm]n  = minimum  of  mass  flow  rate  times  heat  capacity  for  either  fluid,  W/K 
(Btu/(h-°F) 

= maximum  of  the  mass  flow  rate  times  the  heat  capacity  for  either  fluid,  W/K 
(Btu/(h-°F) 

Cpa  = specific  heat  at  constant  pressure  for  air,  kJ/(kg-K)  (Btu/(lb*°F)) 

cfm  = airflow  in  cubic  feet  per  minute,  ft  /min 

COP  = coefficient  of  performance 

fpm  = velocity  in  feet  per  minute,  ft/min 

hj  = inside-tube  heat-transfer  coefficient,  W/(m2-K),  (Btu/(ft2-h-°F)) 

h]  = heat-transfer  coefficient  for  condensate  layer,  W/(m2-K)  (Btu/(ft2*h-°F)) 

hPf  = heat-transfer  coefficient  for  tube/fin  contact,  W/(m2-K)  (Btu/(ft2-h-°F)) 

ho  = air-side  heat  transfer  coefficient,  W/(m2-K)  (Btu/(ft2-h-°F)) 

ifg  = latent  heat  of  evaporation,  kJ/(kg-K)  (Btu/(lb°F) 

in  WG  = inches  of  water  in  gage  pressure 

IP  = inch-pound  or  English  system  of  units 

K = material  thermal  conductivity,  W/(m*K)  (Btu/(ft-h-°F)) 

L = length,  m2  (ft2) 

ma  = air  mass  flow  rate,  kg/s  (lb/h) 

NTU  = number  of  transfer  units  for  the  heat  exchanger,  dimensionless 
Q = capacity  or  heat  transferred,  W (Btu/h) 

scfm  = airflow  in  standard  cubic  feet  per  minute  where  flowrate  is  taken  at  the  air  standard 
density  of  0.075  lbm/ft3  (ANSI/ASHRAE  51-1985) 


SI  = international  system  of  units  or  metric  system  of  units 
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T = temperature,  K (°F) 

t = thickness,  m (ft) 

TXV  = thermostatic  expansion  valve 

U = overall  heat-transfer  coefficient,  W/(m2-K)  (Btu/(ft2-h-°F)) 

W = width,  m (ft) 

Xp  = thickness  of  the  tube  wall,  m (ft) 

OC  — ifgw(tOa  “ t1Xv)/(Cpa(rfa  ~ Tw)) 

£ = heat-transfer  effectiveness,  fraction 

(J)  = fin  efficiency,  fraction 

X = longitudinal  heat  conduction  parameter,  dimensionless 

x = tube  longitudinal  conduction  effect  factor,  fraction 

C0a  = humidity  ratio  of  air  at  tube  inlet,  kgw/kga(dry  (lbw/lba,dry) 

coo  = humidity  ratio  of  saturated  air  at  temperature  of  condensate  wetting  the  tube, 

kgw/kga  ,dry  (lbw/lba>dry) 

Subscripts 
a = air 

f = fin 

i = inlet,  inside,  or  tube  numbering  index 

j = tube  numbering  index 

nc  = no  longitudinal  conduction  effects 

wc  = considering  longitudinal  conduction  effects 

r = refrigerant 

sim  = simulation 

w = tube  wall  or  water 
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1.  SCOPE  OF  THE  STUDY 


Typically,  firmed  tube  evaporators  employ  parallel  refrigerant  circuits  to  obtain  an 
optimal  refrigerant  mass  flux,  which  affects  refrigerant  heat  transfer  coefficient  and 
pressure  drop.  Each  circuit  performs  optimally  when  the  superheat  at  its  exit  matches  the 
desired  overall  superheat  in  the  exit  manifold.  Circuit  superheat  is  affected  by  the 
refrigerant  mass  flowrate  and  the  air  flowrate  associated  with  each  tube. 

Most  evaporators  use  an  inlet  expansion  valve  with  a flow  distributor  to  control  the  bulk 
superheat  at  the  evaporator  exit  manifold.  The  current  practice  does  not  embody  means 
for  adjusting  refrigerant  distribution  between  different  circuits  as  needed.  This  means 
that  non-uniform  airflow  or  unintended  pressure  drops  could  cause  some  circuits  to  have 
excessive  superheat  while  others  may  remain  two-phase  at  the  evaporator  exit.  In  such 
situations,  some  circuits  are  inefficiently  using  coil  area  by  transferring  heat  with 
superheated  vapor  instead  of  two-phase  refrigerant.  There  are  also  mixing  losses 
associated  with  reaching  the  final  superheat  as  two-phase  and  superheated  refrigerant  mix 
in  the  evaporator  exit  manifold. 

The  advances  in  micro  electro-mechanical  systems  (MEMS)  offers  the  opportunity  to 
develop  and  place  inexpensive  flow  control  valves  on  each  circuit  of  an  evaporator.  This 
would  allow  control  of  the  refrigerant  superheat  at  the  exit  of  each  circuit.  This 
experimental  investigation  examines  the  benefits  of  controlling  superheat  by  placing 
individual  needle  valves  on  each  circuit  of  the  evaporators.  The  study  involves  three 
evaporators  containing  three  parallel  refrigerant  circuits  in  identical  configurations.  Two 
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of  these  evaporators  equipped  with  enhanced  (wavy-lanced)  fins,  respectively,  were 
tested  to  examine  the  benefits  of  maintaining  even  superheat  when  compared  to  three 
different  scenarios  of  excessive  superheat.  The  third  evaporator  with  wavy-lanced  fms 
and  separated  (cut)  depth  rows  facilitates  documenting  the  impact  of  tube-to-tube  heat 
conduction.  Non-uniform  superheats  are  imposed  and  compared  to  uniform  superheats  of 
5.6  °C  (10.0  °F)  and  16.7  °C  (30.0  °F).  Non-uniform  airflow  is  also  imposed  while 
superheats  are  allowed  to  adjust  naturally,  and  while  superheats  are  controlled  by  the 
individual  expansion  valves.  The  NIST  tube-by-tube  evaporator  model,  EVAP5,  is  also 
modified  and  used  to  simulate  the  experimental  results. 

The  modeling  part  of  the  study  discusses  longitudinal  tube  and  fin  conduction  and 
presents  a scheme  for  including  tube-to-tube  heat  transfer  into  a tube-by-tube  simulation 
model.  Validated  results  for  an  upgraded  evaporator  model  are  presented. 
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2.  BACKGROUND  AND  LITERATURE  REVIEW 


Refrigerant  incurs  a phase  change  from  the  two-phase  to  the  superheat  zone  in  the 
evaporator.  Large  refrigerant  mass  flux  not  only  increases  the  heat  transfer  rate,  but  also 
increases  the  pressure  drop.  Therefore,  most  refrigerant  evaporators  employ  parallel 
circuits  to  provide  a balanced  effect  on  the  evaporator  capacity  between  the  negative 
effect  of  refrigerant  pressure  drop  and  the  positive  effect  of  improved  inside-tube  heat 
transfer  coefficient. 

Even  though  all  refrigerant  circuits  have  the  same  inlet  and  outlet  conditions,  the 
refrigerant  distribution  is  not  uniform;  the  staggered  tube  arrangement  can  cause  different 
heat  transfer  rates.  Non-uniform  refrigerant  distribution  is  also  due  to  the  thermal 
resistance  in  the  superheat  region  increasing  more  rapidly  than  in  the  two-phase  region. 
Superheated  vapor  at  the  evaporator  exit  is  necessary  to  prevent  liquid  compression  and 
subsequent  damage  to  the  compressor,  even  though  superheat  reduces  the  performance  of 
the  system. 

Refrigerant  superheat  in  a given  circuit  is  affected  by  the  refrigerant  mass  flow  rate  and 
the  airflow  rate  over  the  coil  area  associated  with  that  circuit.  For  a given  air  distribution 
there  is  one  refrigerant  flow  rate  that  results  in  a desired  superheat  at  the  individual 
circuit  exit.  When  circuits  are  not  well  balanced,  the  target  overall  superheat  is  a result  of 
mixing  a highly  superheated  refrigerant  and  two-phase  refrigerant  leaving  different 
circuits.  This  causes  significant  degradation  in  evaporator  capacity  because  the  circuit 
with  superheated  refrigerant  transfers  less  heat. 
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Liang  et  al.(2001)  conducted  a numerical  study  of  the  refrigerant  circuit.  The  governing 
equations  and  control  volumes  were  presented  with  the  simulation  procedure  for 
branches,  tubes,  and  control  volumes  of  a coil.  Using  the  model,  the  heat  transfer  and 
fluid  flow  characteristics  of  the  coils  were  studied.  Compared  to  a common  coil,  the 
researchers  found  that  using  a complex  refrigerant  circuit  arrangement  where  the 
refrigerant  circuits  are  properly  branched  or  joined  may  reduce  the  heat  transfer  area  by 
around  5 % while  maintaining  constant  capacity.  The  investigators  experimentally 
validated  6 different  refrigerant  circuiting  arrangements  while  maintaining  the 
evaporators  inlet  and  exit  states.  They  used  an  R134a  expansion  valve  inlet  condition  of 
40  °C  (104  °F)  saturated  liquid  with  5.0  °C  (2.8  °F)  of  subcooling  and  an  evaporator  exit 

saturation  temperature  of  10  °C  (50  °F)  with  5.0  °C  (2.8  °F)  of  superheat.  Liang  et  al. 
noted  that  for  a given  evaporator  load,  designers  must  design  the  refrigerant  circuitry  to 
produce  a refrigerant  mass  velocity  that  produces  a maximum  heat  flux.  Maximum  heat 
fluxes  vary  with  refrigerant  circuiting  due  to  varying  levels  of  refrigerant  pressure  drop. 
Their  model  was  able  to  predict  evaporator  capacity  within  5 % on  four  of  the  six  coils 
while  predicting  refrigerant  pressure  drop  to  within  25  %. 

Kirby  et  al.  (1998)  experimentally  investigated  the  performance  of  a 5275  W 
(18  000  Btu/h)  window  air  conditioner  under  wet  and  dry  coil  conditions  with  non- 
uniform  airflow  over  the  evaporator.  The  velocity  variation  over  the  evaporator  varied  by 
as  much  as  a factor  of  3,  but  upon  correcting  the  non-uniformity  of  airflow,  the 
investigators  saw  only  a minor  improvement  in  performance.  This  was  a system  study 
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with  no  attempt  to  maintain  constant  refrigerant  states  at  the  inlet  and  exit  of  the 
evaporator.  A round  disk  was  used  to  block  16  % of  the  central  area  of  the  evaporator 
while  maintaining  the  original  airflow.  Tests  were  conducted  with  this  blockage  against 
the  evaporator  face  and  then  moved  in  steps  in  the  upstream  direction.  They  found  no 
capacity  degradation  greater  than  2 %.  The  authors  noted  that  the  blockage  caused  more 
of  the  evaporating  refrigerant  to  exist  in  the  two-phase  state;  thereby  reducing  the 
superheated  area  of  the  evaporator  and  offsetting  the  inability  of  the  air  velocity  to 
compensate  for  the  loss  of  heat  transfer  area.  Wet-coil  tests  also  showed  very  little 
difference  in  the  sensible  heat  ratio  with  non-uniform  airflow.  The  authors  noted  that  any 
non-uniformity  in  airflow  must  be  noted  by  designers  and  used  to  intelligently  circuit  the 
refrigerant  to  equalize  exposure  of  the  evaporating  refrigerant  to  airflow. 

Chwalowski  et  al.  (1989)  examined  computer  models  and  performed  experiments  using 
three  different  evaporators;  two  V-shaped  evaporators  with  upflow  and  one  vertical  slab 
evaporator  with  horizontal  and  angled  flow  with  respect  to  the  approaching  airstream. 
These  were  evaporator  tests  with  fixed  evaporator  saturation  pressures;  therefore,  these 
tests  parallel  the  technique  used  in  the  current  experimental  investigation.  The 
investigators  determined  coil  face  velocity  for  several  of  the  configurations  and  noted 
non-uniformity  of  the  airflow.  Generally,  the  evaporator  capacity  varied  with  the 
configuration  mainly  as  a function  of  the  exit  superheats  at  the  two  evaporator  exits.  As 
the  coil  angle  with  respect  to  the  approaching  airstream  was  varied,  capacity  degradations 
on  the  order  of  20  % were  noted.  Non-uniformities  in  superheat  were  produced  by  non- 
uniform  airflow  that  produced  differences  in  heat  transfer  and  pressure  drops  within  the 
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heat  exchangers.  The  design  of  the  circuitry  and  various  splitting  points  within  the 
evaporator  produced  differing  capacities  as  a function  of  the  coil  orientation  with  the 
airstream.  The  investigators  noted  that  none  of  the  three  evaporator  models  they 
considered  could  accurately  predict  capacity  without  some  knowledge  of  the  air  velocity 
profile  and  refrigerant  maldistribution. 
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3.  LABORATORY  EXPERIMENT 


3.1  Experimental  Setup 

Figure  3.1.1  shows  a schematic  diagram  of  the  experimental  setup.  The  test  rig  consisted 
of  three  major  flow  loops:  (1)  a refrigerant  flow  loop  containing  a detachable  test 
section,  (2)  a water  flow  loop  used  for  the  condensation  heat  exchanger  and  (3)  an  air 
flow  loop  used  for  the  evaporation  heat  exchanger.  The  design  of  the  rig  allowed  easy 
control  of  operating  parameters  such  as  condensing  pressure  and  subcooling  at  the  inlet 
of  the  expansion  valve  (evaporator  inlet  enthalpy),  evaporating  pressure  at  the  exit  of  the 
evaporator,  and  superheat  at  the  exit  of  the  evaporator. 


(2)  Water  Flow  Loop 
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Figure  3.1.1 : Schematic  diagram  of  the  experimental  setup 
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Figure  3.1.2  is  a photo  of  the  specially  designed  and  constructed  R22  condensing  unit. 
The  design  of  this  condensing  unit  allowed  complete  control  of  the  subcooled  R22  liquid 
conditions. 


Figure  3.1.2:  Condensing  unit  used  to  precisely  control  refrigerant  conditions 


An  open-drive  compressor  with  a variable  speed  motor  provided  refrigerant  mass  flow, 
set  the  enthalpy  entering  the  test  section,  set  the  condensing  pressure,  and  set  the 
subcooling  at  the  inlet  of  the  expansion  valves.  We  controlled  condensing  pressure  by 
adjusting  cooling  water  flow  using  a hand-operated  needle  valve.  To  provide  additional 
pressure  control  for  the  condenser,  we  also  controlled  the  entering  temperature  of  the 
cooling  water.  Water  flow  rate  and  temperature  through  the  subcooler  plate  type  heat 
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exchanger  controlled  the  refrigerant  subcooling  at  the  inlet  of  the  expansion  valve.  A 
loop  supplying  water  to  the  condensing  heat  exchanger  and  subcooler  consisted  of  a 
refrigerator,  water  storage  tank,  and  a pump.  The  controls  combined  to  produce 
evaporator  inlet  quality  of  25  % ± 1 % (using  a 95  % confidence  interval). 

The  test  section  has  three  parallel  refrigerant  paths.  The  exit  pressure  of  the  evaporator 
was  adjusted  by  a pressure  regulating  valve  which  was  installed  in  the  refrigerant  line. 
The  superheats  of  the  three  circuits  of  the  evaporator  were  adjusted  by  three  manual 
expansion  valves.  The  design  of  the  test  section  allowed  easy  installation  and 
replacement  of  the  evaporators.  The  pressure  and  temperature  were  measured  upstream 
and  downstream  of  the  test  rig.  Flow  conditions  were  also  monitored  using  a sight  glass 
at  the  inlet  of  the  compressor  and  expansion  valves. 

The  total  refrigerant  flow  rate  was  measured  by  a coriolis-type  mass  flow  meter  in  the 
liquid  line  between  the  subcooler  and  the  expansion  valves.  Two  turbine  flow  meters 
were  installed  to  measure  flow  rate  in  two  of  the  circuits.  Flow  through  the  third  circuit 
was  calculated  by  subtracting  the  flow  through  two  of  the  circuits  from  the  total  mass 
flow. 

Air  flow  rate  was  measured  in  the  air  flow  chamber  according  to  ANSI/ASHRAE  51- 
1985.  Evaporator  capacity  was  calculated  using  the  air  enthalpy  method  and  refrigerant 
enthalpy  method  following  procedures  specified  in  ASHRAE  Standard  37  (1998).  In  the 
present  experiments,  the  maximum  difference  between  the  air  and  refrigerant  side 
capacity  was  less  than  5 %.  Air  velocity  at  the  face  of  the  evaporator  was  measured  with 
a hot  wire  anemometer. 
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Table  3.1.1  lists  the  parameters  controlled  to  produce  a successful  evaporator  test.  The 
parameters  are  listed  in  the  order  they  were  set  to  produce  a controlled  test. 

Table  3.1.1:  Essential  Control  Parameters 


Parameter 

Setpoint 

Upstream  Liquid  Saturation  Temperature 

40.6  °C  (105.0  °F):  controlled  by  condenser 
water  flowrate  and  compressor  speed 

Liquid  Line  Subcooling 

8.3  °C  (15.0  °F):  controlled  by  upstream 
pressure  and  refrigerant  charge 

Evaporator  Circuit  Superheats 

5.6  °C  or  16.7  °C  (10.0  °F  or  30.0  °F): 
controlled  by  expansion/needle  valve 
opening  and  evaporator  exit  pressure 

Evaporator  Exit  Saturation  Temperature 

7.2  °C  (45.0  °F):  controlled  by  evaporator 
pressure  regulator  valve  and  compressor 
speed 

Evaporator  Inlet  Liquid  Enthalpy 

Corresponds  to  the  saturated  liquid 
temperature  of  40.6  °C  (105.0  °F):  when 
inlet  pressures  were  increased,  the  inlet 
enthalpy  was  always  monitored  to  produce 
an  enthalpy  equal  to  the  saturated  liquid 
enthalpy  at  40.6  °C  (105.0  °F)  ± 1.4  °C 
(2.5  °F) 
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3.2  Evaporators  Selected  for  Testing 

We  used  three  finned  tube  heat  exchangers  of  the  same  outside  dimensions,  tube  spacing, 
and  circuitry  as  the  test  evaporators:  (1)  COIL-W  with  wavy  fins,  (2)  COIL-E  with  wavy- 
lanced  (enhanced)  fins,  and  (3)  COIL-EC  (Figure  3.2.1)  with  wavy-lanced  (enhanced) 
fins  and  the  tube  rows  separated  to  inhibit  tube-to-tube  heat  transfer  (enhanced-cut). 
Figures  3.2.2,  3.2.3,  and  3.2.4  show  the  side  views  of  the  refrigerant  circuits.  The 
following  are  the  main  design  parameters: 

(a)  3 depth  rows  with  25.4  mm  (1  in)  face  spacing  and  22.0  mm  (0.866  in)  row 
spacing 

(b)  3 refrigerant  circuits  as  shown  in  Figures  3.2.1  and  3.2.2 

(c)  9.53  mm  (0.375  in)  diameter  round  copper  tubes,  smooth  walls,  0.254  mm 
(0.010  in)  wall  thickness 

(d)  0.1 143  mm  (0.0045  in)  thick  aluminum  fins;  wavy  fins  for  COIL-W  and 
louvered  or  slit  fins  for  COIL-E  and  COIL-EC 
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Figure  3.2.1:  COIL-EC  showing  separated  tube  depth  rows 
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Figure  3.2.2:  A schematic  side  view  of  refrigerant  circuitry 
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Figure  3.2.3:  A schematic  side  view  of  refrigerant  circuitry  for  COIL-EC 


Inlets 


Figure  3.2.4:  Circuiting  of  all  three  evaporators 
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3.3  Test  Conditions  and  Experimental  Procedure 

Table  3.3.1  lists  the  test  parameters  and  environmental  chamber  conditions  for  tests  on 
the  three  evaporators.  All  tests  were  conducted  at  the  same  26.7  °C  (80.0  °F)  indoor  dry- 

bulb  with  dew-point  varying  for  wet-coil  tests  and  dry-coil  tests.  Refrigerant  R22 
conditions  at  the  inlet  to  the  expansion  valves  were  controlled  to  maintain  an  enthalpy 
equivalent  to  a 48.9  °C  (120.0  °F)  saturation  temperature  with  a subcooling  of  8.3  °C 

(15.0  °F)  ± 1.4  °C  (2.5  °F).  Some  tests  required  increasing  the  inlet  pressure  to  produce 
the  required  superheats  at  the  evaporator  circuit  exits.  When  the  pressure  was  increased, 
the  enthalpy  and  subcooling  were  adjusted  to  keep  a constant  enthalpy  at  the  evaporator 
inlet. 


Table  3.3.1 : Experimental  Test  Conditions 


Variable 

Value 

Tolerance 

Indoor  Dry- Bulb 

26.7  °C  (80.0  °F) 

0.28  °C  (0.5  °F) 

Indoor  Dew-Point  for  Wet- 
Coil  Tests 

15.8  °C  (60.4  °F) 

0.28  °C  (0.5  °F) 

Evaporator  Exit  Saturation 
Temperature 

7.2  °C  (45.0  °F) 

0.28  °C  (0.5  °F) 

Evaporator/Expansion 
Valve  Inlet  Saturation 
Temperature 

48.9  °C  (120.0  °F) 

1.4  °C  (2.5  °F) 

Evaporator/Expansion 
Valve  Inlet  Subcooling 

8.3  °C  (15.0  °F) 

1.4  °C  (2.5  °F) 

Once  an  evaporator  coil  was  mounted  in  the  airflow  chamber,  indoor  dry-bulb  and  dew- 
point were  stabilized  for  at  least  1 h.  While  indoor  psychrometric  conditions  stabilized, 
the  evaporator  inlet  R22  pressure  and  temperature  were  set  by  adjusting  the  flow  control 
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valves  on  the  condensing  unit.  Water  flow  to  the  condenser  and  subcooler  plate  heat 
exchangers  was  adjusted  to  establish  the  evaporator  expansion  valves  inlet  pressure  and 
temperature.  The  evaporator  exit  saturation  temperature  was  set  by  adjusting  the 
evaporator  pressure  regulating  valve  at  the  exit  of  the  evaporator.  Superheat  conditions 
in  the  individual  circuits  were  set  by  adjusting  R22  mass  flow  through  each  circuit. 
Airflow  rate  over  the  evaporator  was  adjusted  using  the  variable  speed  drive  on  the 
airflow  chamber’s  pull-thru  fan. 

Table  3.3.2  and  3.3.3  list  the  tests  performed  for  the  evaporators.  Capacity  specific 
airflow  rate  was  initially  established  at  193  nT/kW  h (400  scfm/ton)  for  test  9.  Test  9 
required  manipulating  superheats  and  airflow  rate  to  obtain  the  desired  airflow  to 
capacity  ratio.  Test  9 capacity  was  then  used  to  calculate  the  airflow  rate  for  the 
169  m3/kW  h (300  scfm/ton)  and  242  mJ/kW  h (500  scfm/ton)  tests.  Tests  with  a cross- 
parallel flow  configuration  were  performed  by  switching  the  refrigerant  flow. 

Non-uniform  airflow  tests  were  performed  with  COIL-W  and  COIL-E.  We  established 
non-uniform  air  distribution  by  attaching  a series  of  metal  mesh  plates  to  the  upper  half 
of  the  coil. 

A hot  wire  anemometer  was  used  to  measure  airflow  rate  by  traversing  the  coil  at  a 
minimum  of  25  equally  spaced  points  at  the  face  of  the  coil.  This  measurement  agreed 
with  the  chamber  airflow  within  2 %. 
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Table  3.3.2:  Evaporator  Tests 


Evaporator 

Tests  Performed 

COIL-W  (wavy  fins),  cross-counter  flow 

1,  5-13  (10  tests) 

COIL-W  (wavy  fins),  cross-parallel  flow 

9-12  (4  tests) 

COIL-E  (lanced  fins),  cross-counter  flow 

1,2,  5-14  (12  tests) 

COIL-E  (lanced  fins),  cross-parallel  flow 

9-12  (4  tests) 

COIL-EC,  cross-counter  flow 

1,2,  5,  6,  9,  10,  13,  14  (8  tests) 

COIL-W,  cross-counter  flow,  non-uniform 
airflow 

9 (1/2  profile,  no  superheat  adjustment),  9(1/2 
profile,  superheat  adjusted),  9(1/3  profile,  no 
superheat  adjustment),  9 (1/3  profile,  superheat 
adjusted)  [4  tests] 

COIL-E,  cross-counter  flow,  non-uniform 
airflow 

9(1/2  profile,  no  superheat  adjustment),  9(1/2 
profile,  superheat  adjusted),  9(1/3  profile,  no 
superheat  adjustment),  9(1/3  profile,  superheat 
adjusted)  [4  tests] 

Table  3.3.3:  Test  Number  and  Conditions  for  Each  Evaporator  Test 


% 

vs 

a 

H 

Volumetric  Flowrate 
of  Air  m3/h  (scfm) 

Coil  Surface 

Overall  Superheat 

Superheats  in  Individual  Circuits 

145-Q1 

(300Q) 

193-Q1 

(400-Q) 

242-Q1 

(500-Q) 

>> 

Q 

4-* 

0-> 

£ 

5.6  °C 
(10.0  °F) 

16.7  °C 
(30.0  °F) 

5.6  °C 
(10.0  °F) 

5.6  °C 
(10.0  °F) 

5. 6/5. 6/5. 6 
(10/10/10) 

16.7/16.7/16. 

7 (30/30/30) 

16.7/*/16.7 

(30/*/30) 

*/16.7/16.7 

(*/30/30) 

1 

X 

X 

1 

2 

X 

X 

2 

3 

X 

X 

3 

4 

X 

X 

4 

5 

X 

X 

5 

6 

X 

X 

6 

7 

X 

X 

7 

8 

X 

X 

8 

9 

X 

X 

9 

10 

X 

X 

10 

11 

X 

X 

11 

12 

X 

X 

12 

13 

X 

X 

13 

14 

X 

X 

14 

15 

X 

X 

15 

16 

X 

X 

16 

* Superheat  to  be  controlled  such  that  the  desired  overall  level  of  superheat  is  obtained 
1)  SI  units  of  m3/kW  h multiplied  by  capacity  (Q)  in  kW  to  determine  airflow,  (IP  units  of  cfm/ton 
multiplied  by  capacity  (Q)  in  tons). 
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In  total  we  performed  54  tests  with  uniform  airflow  and  28  tests  with  imposed  non- 
uniform  distribution  of  air.  We  conducted  a total  of  90  tests  including  repeats  and  tests 
that  were  excluded  due  to  unsteady  or  non-standard  conditions. 

The  capacity  characteristics  of  the  three  evaporators  are  shown  below.  Figure  3.3.1 
shows  the  capacity  ratio  at  different  test  conditions  to  the  capacity  at  test  9 for  the  wavy 
coil.  Tests  1,  9,  and  13  are  wet  coil  tests,  and  test  5 is  a dry  coil  test.  COIL-E  and  COIL- 
EC  evaporators  represented  higher  capacity  than  that  of  the  COIL-W  evaporator  for  the 
wet  coil  tests.  Even  though  the  air-side  sensible  heat  transfer  coefficient  is  much  lower 
than  the  refrigerant  side,  the  air-side  thermal  resistance  is  reduced  due  to  the 
enhancement  of  moisture  condensation  and  large  finned  area.  For  wet  coil  tests  (tests  1, 
9,  13),  the  capacity  of  COIL-E  and  COIL-EC  was  larger  than  that  of  COIL-W.  The 
COIL-EC  produced  higher  capacity  than  COIL-E  possibly  because  of  the  added  fin 
leading  edges  agitating  the  boundary  layer  and  increasing  the  air-side  heat  transfer 
coefficient  even  further  in  addition  to  eliminating  some  tube-to-tube  heat  transfer. 
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Capacity  Ratio  with  respect  to  Wavy  T est  9 


1.40 
1.20 
1.00 
0.80 
0.60 
0.40 
0.20 
0.00 

Figure  3.3.1 : Capacity  ratio  for  different  shape  fins  relative  to  the  capacity  at  test  9 for 
the  wavy  coil  (Test  1 : low  airflow, wet  coil,  Test  5:  median  airflow,  dry 
coil,  Test  9:  median  airflow,  wet  coil,  Test  13:  high  airflow,  wet  coil). 
All  tests  have  a 5.6  °C  (10.0  °F)  uniform  superheat. 
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□ Wavy 

■ Enhanced 

□ Enhanced-Cut 
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3.4  Experimental  Results 

3.4. 1 Cross-Counter  Air/Refrigerant  Flow  Configuration  Tests 

3.4.J.1  Non-Uniform  Superheat  Tests 

Figure  3. 4. 1.1.1  shows  capacity  at  different  superheat  test  conditions.  These  data  are 
shown  in  Table  3.4. 1.1.1.  All  of  the  coils  showed  a rapid  decrease  in  capacity  when 
individual  circuit  superheat  was  increased  with  the  overall  superheat  maintained  at  5.6  °C 
(10  °F).  Figure  3.4. 1.1. 2 shows  the  relative  capacity  of  COIL-W  tests  10  and  12  with 
respect  to  test  9.  Even  though  the  overall  superheat  was  held  at  5.6  °C  (10  °F),  the  non- 
uniformity of  superheat  in  cases  1 1 and  12  produced  a 41  % loss  in  capacity.  This  was 
almost  as  severe  as  the  43  % loss  in  capacity  seen  when  overall  superheat  was  held  at 
16.7  °C  (30.0  °F). 

Test  12  showed  similar  capacity  as  test  10  even  though  overall  superheat  was  lower.  At 
test  condition  12,  the  mass  flow  rate  through  the  top  circuit  was  much  higher  than  that  at 
test  10.  Therefore,  the  inlet  refrigerant  temperature  of  the  evaporator  for  test  10  was 
higher  than  test  12  because  exit  pressure  was  the  same.  This  means  that  the  temperature 
difference  between  air  and  refrigerant  for  test  12  was  higher  than  for  test  10.  This 
allowed  test  12  to  have  a higher  capacity  than  test  10. 
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Table  3. 4. 1.1.1:  Non-Uniform  Superheat  Test  Data  for  COIL-W  and  COIL-E 


Test  Name 

Test  # 

1 

Coil  Designation  jj 

Volumetric  Flowrate  of 
Air  m3/h  (cfm) 

Coil 

Surface 

Overall  Superheat 

Superheats  in  Individual  Circuits 

145Qla 

(300-Q) 

193Qla 

(400-Q) 

242  Qla 
(500-Q) 

Dry 

Wet 

5.6  °C 
(10.0  °F) 

16.7  °C 
(30.0  °F) 

5.6  °C 
(10.0  °F) 

5.6  °C(10.0  °F) 

5. 6/5. 6/5. 6 
(10/10/10) 

16.7/16.7/16.7 

(30/30/30) 

16.7/*/16.7 

(30/*/30) 

*/l  6.7/1 6.7 
(*/30/30) 

Qtest  ^ 

^tes/ 

3lb 

Qtest 

Qtes/ 

Qlb 

Qtest  | 

3test 

Qtest 

Qtest/ 

Qlb 

W (Btu/h) 

W (Btu/h) 

W (Btu/h) 

'Qlb' 

V (Btu/h) 

Cross-Counter  Air/Refrigerant  Flow 

W020225B 

5 

W 

X 

X 

5428 

(18519) 

1 

W020228A 

6 

W 

X 

X 

3569 

(12177) 

0.66 

W020221 A 

7 

w 

X 

X 

3910 

(13341) 

0.72 

W020225A 

8 

w 

X 

X 

3888 

(13266) 

0.72 

W020207B 

9 

w 

X 

X 

6507 

(22203) 

1 

W020530A 

10 

w 

X 

X 

3722 

(12701) 

0.57 

W020531 A 

11 

w 

X 

X 

3837 

(13091) 

0.59 

W020215B 

12 

w 

X 

X 

3830 

(13067) 

0.59 

W020322A 

5 

E 

X 

X 

5602 

(19115) 

1 

W020321B 

6 

E 

X 

X 

4301 

(14677) 

0.77 

W020322C 

7 

E 

X 

X 

4797 

(16367) 

0.86 

W020322B 

8 

E 

X 

X 

4700 

(16037) 

0.84 

E020607A 

9 

E 

X 

X 

6955 

(23733) 

1 

W020318A 

10 

E 

X 

X 

4865 

(16599) 

0.70 

W020318B 

11 

E 

X 

X 

5485 

(18715) 

0.79 

W020319A 

12 

E 

X 

X 

4735 

(16157) 

0.68 

* Superheat  to  be  controlled  such  that  the  desired  overall 


evel  of  superheat  is  obtained 


la)  SI  units  of  m /kW  h multiplied  by  capacity  (Q)  in  kW  to  determine  airflow,  (IP  units  of  cfm/ton 


multiplied  by  capacity  (Q)  in  tons). 


lb)  Capacity  relative  to  the  5.6  °C  (10.0  °F)  tests  noted  by  a ratio  of  1 in  the  row  above. 
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Test  Number 

Figure  3 .4. 1 . 1 . 1 : Capacity  of  evaporators  for  tests  5 through  1 2 


1.2 


Test  9 Test  10  Test  12 

Figure  3.4. 1.1. 2:  Capacity  ratio  at  different  superheat  conditions  relative  to  test  9 for 
COIL-W  and  COIL-E  (superheat  cases  follow  Table  3.4. 1.1.1). 
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Evaporator  Capacity  (W) 


For  COIL-E,  the  reduction  in  capacity  due  to  an  increase  in  superheat  was  lower  than 
COIL-W  (Figure  3 .4. 1.1.1  and  3. 4. 1.1. 2).  COIL-E  seemed  to  show  a preference  as  to 
which  flooded  circuit  (test  11  or  12)  produced  the  higher  capacity.  When  the  middle 
circuit  was  flooded,  the  capacity  decreased  by  21  % compared  to  32  % when  the  top 
circuit  was  flooded. 

Figure  3. 4. 1.1. 3 shows  the  relative  capacity  of  COIL-W  for  tests  10  and  12  with  respect 
to  test  9 and  for  tests  6 and  8 with  respect  to  test  5.  The  capacity  of  the  dry  coil 
decreased  by  28  % with  non-uniform  circuit  superheats  with  overall  superheat  fixed  at 
5.6  °C  (10  °F).  COIL-E  capacity  (Figure  3. 4. 1.1. 4)  dropped  by  16  % with  non-uniform 
superheat  under  dry  conditions;  again  COIL-E  showed  that  flooding  the  middle  circuit 
produced  a smaller  capacity  drop  than  flooding  the  top  circuit  while  holding  overall 
superheat  constant  at  5.6  °C  (10  °F). 
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Test  9 and  5 Test  10  and  6 Test  12  and  8 


Figure  3.4. 1 . 1 .3:  Capacity  ratio  for  COIL-W  relative  to  test  9 for  wet  and  test  5 for  dry 
conditions 
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Test  9 and  5 Test  10  and  6 Test  12  and  8 


Figure  3.4. 1 . 1 .4:  Capacity  ratio  for  COIL-E  relative  to  test  9 for  wet  and  test  5 for  dry 
conditions 


The  experimental  investigation  was  designed  to  reveal  some  of  the  effects  of  tube-to-tube 
heat  transfer  by  heat  conduction  through  the  fin  material.  The  comparison  was  done  by 
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examining  COIL-E  and  COIL-EC  and  comparing  their  capacity  at  different  levels  of 
superheat.  In  addition  to  the  capacity  comparison  between  COIL-E  and  COIL-EC,  direct 
evidence  of  conduction  between  tubes  was  noted  from  the  thermocouple  bend 
temperature  data  for  COIL-W. 

Figures  3.4. 1.1.5  and  3.4. 1.1. 6 show  bend  temperature  data  for  test  9 (6.7  °C  (10  °F) 
superheat  on  all  circuits)  and  test  12  (flooded  top  circuit  with  16.7  °C  (30.0  °F)  superheat 

on  the  other  two  circuits).  The  figure  for  test  9 shows  a uniform  temperature  distribution 
between  comparable  bends  in  the  three  refrigerant  circuits.  The  noticeable  difference 
occurs  for  test  12  where  the  low  temperature,  flooded  top  circuit  shows  some  thermal 
communication  with  the  final  tube  passes  of  the  middle  circuit.  The  top  circuit  is 
showing  an  average  surface  temperature  of  approximately  9.4  °C  (49  °F)  throughout  all 
of  its  tubes,  while  the  middle  circuit  shows  definite  superheat  at  the  third  and  fourth  final 
tube  bend  with  a temperature  of  24.0  °C  (75.2  °F).  This  is  where  the  conduction  between 

circuits  was  obvious;  the  surface  temperature  on  the  final  two  tubes  bend  was  22.3  °C 
(72.2  °F).  This  is  a decrease  in  temperature  due  to  conduction  between  the  top  circuit’s 
tube  and  the  middle  circuit’s  tubes. 

The  conduction  effects  were  quantified  in  the  tests  conducted  with  COIL-E  and  COIL- 
EC;  by  separating  the  tube  sheets  in  COIL-EC  and  thereby  removing  a majority  of  the 
conduction  path  between  tubes.  Figure  3.4. 1.1. 7 shows  the  capacity  of  COIL-E  and 
COIL-EC  relative  to  test  9 for  COIL-E  during  cross-counter  flow  for  tests  9,  10,  13,  and 
14.  These  two  coils  used  identical  fin  material  and  fin  type;  the  only  difference  was  the 
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tube  sheets  of  COIL-EC  were  separated.  Figure  3.4. 1.1. 8 shows  that  for  test  9,  COIL-E 
and  COIL-EC  have  very  similar  bend  temperatures.  Figure  3. 4. 1.1. 9 shows  the  same 
coils  with  the  superheat  increased  to  16.7  °C  (30.0  °F)  at  the  coil  exit  (test  10).  COIL-EC 
shows  lower  inlet  temperatures  than  COIL-E  even  though  expansion  valve  inlet  and  coil 
exit  conditions  are  almost  identical  for  both  tests.  The  differences  in  temperatures  seen 
with  tests  9 and  10  are  more  pronounced  for  tests  13  and  14  at  the  higher  airflow  rate. 

These  differences  in  temperatures  could  have  produced  the  differences  in  capacity  seen 
between  COIL-E  and  COIL-EC.  As  noted  above,  the  inlet  and  exit  conditions  for  these 
coils  were  nearly  identical,  but  COIL-EC  always  showed  lower  bend  temperatures  than 
COIL-E.  This  would  mean  that  COIL-EC  was  operating  at  a higher  average  temperature 
difference  with  respect  to  the  air  than  COIL-E.  The  greater  average  temperature 
difference  for  COIL-EC  could  translate  to  higher  capacity  than  COIL-E.  The  test  results 
showed  that  both  coils  produced  very  similar  capacities  when  the  overall  superheat  was  at 
5.6  °C  (10.0  °F).  As  the  superheat  was  increased  to  16.7  °C  (30.0  °F),  COIL-EC  capacity 
was  10  % higher  than  COIL-E.  As  the  airflow  increased  for  tests  13  and  14,  COIL-E  and 
COIL-EC  still  produced  nearly  equal  capacities  at  5.6  °C  (10.0  °F)  superheat,  but  when 
superheat  was  increased  to  16.7  °C  (30.0  °F),  COIL-EC  had  a 23  % higher  capacity  than 
COIL-E  (Figures  3.4.1.1.10  and  3.4.1.1.11).  This  tends  to  lend  more  evidence  to 
conduction  effects  between  the  tube  sheets;  eliminating  some  conduction  paths  improved 
the  performance  of  the  enhanced  fin  coil. 
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COIL-W,  Test  9.  Bend  Temperature  (°F) 


Exit  Refrigerant 
^ 54  6 


Figure  3.4. 1 . 1 .5:  COIL-W,  test  9,  circuit  bend  temperatures  for  cross-counter  flow  (All 
circuits  with  6.7  °C  (10.0  °F)  superheat;  refrigerant  exit  manifold 
saturation  temperature  set  to  7.2  °C  (45.0  °F)) 


COIL-W.  Test  12,  Bend  Temperatures0^) 


Exit  Refrigerant 
> 47  8 


► 72.2 


► 73  8 


> 

o' 

$ 


Figure  3.4. 1.1. 6:  COIL-W,  test  12,  circuit  bend  temperatures  for  cross-counter  flow  (Top 
circuit  flooding  with  16.7  °C  (30.0  °F)  superheat  on  bottom  two  circuits 
to  yield  overall  exit  superheat  of  5.6  °C  (10.0  °F);  refrigerant  exit 
manifold  saturation  temperature  set  to  7.2  °C  (45.0  °F)) 
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Figure  3.4. 1.1. 7:  Capacity  of  COIL-E  and  COIL-EC  relative  to  COIL-E,  test  9 at  two 
different  airflow  rates 


Exit  Refrigerant 
► 55.8 


► 56.0 


► 55.8 


Figure  3.4. 1.1. 8:  COIL-EC  and  COIL-E  bend  temperatures  for  test  9 (cross-counter 
flow,  wet  coil,  5.6  °C  (10.0  °F)  superheat  on  all  circuits) 
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Airflow 


Figure  3. 4. 1.1. 9:  COIL-EC  and  COIL-E  bend  temperatures  for  test  10  (cross-counter 
flow,  wet  coil,  16.7  °C  (30.0  °F)  superheat  on  all  circuits) 


Figure  3.4.1.1.10:  COIL-EC  and  COIL-E  bend  temperatures  for  test  13  (cross-counter 

flow,  wet  coil,  5.6  °C  (10.0  °F)  superheat  on  all  circuits) 
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Airflow 


Figure  3.4.1.1.11:  COIL-EC  and  COIL-E  bend  temperatures  for  test  1 4 (cross-counter 

flow,  wet  coil,  16.7  °C  (30.0  °F)  superheat  on  all  circuits) 


3. 4. 1.2  Effects  of  Airflow  Rate  on  Coil  Capacity 

Air  flow  rate  through  the  evaporator  plays  an  important  role  in  the  capacity.  When  air 
flow  rate  is  higher  than  the  optimum  quantity,  the  COP  of  the  system  decreases  due  to  an 
increase  in  the  air  pressure  drop  and  accompanying  fan  power  consumption.  But  higher 
airflow  enhances  the  heat  transfer  rate  of  the  evaporator  on  the  air-side  due  to  the  higher 
Reynolds  number.  Table  3.4. 1.2.1  and  figure  3. 4. 1.2.1  show  the  capacity  variation  of  the 
tested  evaporators  as  a function  of  air  flow  rate.  As  the  air  flow  rate  increased,  total 
capacity  increased  due  to  the  increased  air  mass  flow  rate.  The  latent  heat  transfer  rate 
changed  a little  due  to  the  constant  evaporator  pressure  set  by  the  evaporator  pressure 
regulating  valve.  Other  reasons  that  could  play  a small  part  in  the  constant  latent 
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capacity  may  be  the  condensed  water  on  the  surface  of  tube  and  fin  mixing  with  the  air 
which  is  unsaturated  before  latent  heat  transfer  takes  place  at  the  range  of  these  air  flow 
rates.  Secondly,  the  temperature  difference  between  the  air  and  the  surface  of 
condensing  water  decreases  because  thermal  resistance  increases  due  to  the  condensed 
water  layer.  As  a result,  the  dominant  increase  in  total  capacity  was  caused  by  the 
increase  in  the  sensible  heat  transfer  rate. 


Table  3.4. 1.2.1:  Capacity  of  the  Test  Evaporators  at  Varying  Airflow  Rates 


Test 

Name 

Coil  Designation  || 

Test  # jj 

Volumetric  Flowrate 
of  Air  m3/h  (scfm) 

Coil  Surface 

Overall  Superheat 

Superheats  in  Individual  Circuits 

145-Q1 

(300-Q) 

193-Q1 

(400-Q) 

242-Q1 

(500-Q) 

>> 

i— 

Q 

<u 

£ 

5.6  °C 
(10.0  °F) 

16.7  °C 
(30.0  °F) 

5.6  °C 
(10.0  °F) 

5.6  °C 
(10.0  °F) 

5. 6/5. 6/5. 6 
(10/10/10) 

16.7/16.7/16.7 

(30/30/30) 

16.7/*/16.7 

(30/*/30) 

*/l  6.7/1 6.7 
(*/30/30) 

W (Btu/h) 

W (Btu/h) 

W (Btu/h) 

W (Btu/h) 

W020226A 

W 

1 

X 

X 

5788 

(19746) 

W020207B 

W 

9 

X 

X 

6508 

(22203) 

W020301 A 

w 

13 

X 

X 

7503 

(25598) 

W020320B 

E 

1 

X 

X 

5998 

(20464) 

W020607A 

E 

9 

X 

X 

6956 

(23732) 

W020319B 

E 

13 

X 

X 

7653 

(26109) 

E020417A 

EC 

1 

X 

X 

6085 

(20760) 

E020415A 

EC 

9 

X 

X 

6972 

(23788) 

E020416B 

EC 

13 

X 

X 

7781 

(26546) 

* Superheat  to  be  controlled  such  that  the  desirec 


overall  level  of  superheat  is  obtained 


1)  SI  units  of  m /kWh  multiplied  by  capacity  (Q)  in  kW  to  determine  airflow,  (IP  units  of  cfm/ton 


multiplied  by  capacity  (Q)  in  tons). 
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Air  flow  rate(CFM) 


Figure  3.4. 1.2.1 : Capacity  as  a function  of  air  flow  rate  for  wet  coil  conditions 


3. 4. 1.3  Effects  of  Non-Uniform  Airflow  on  Coil  Capacity 

The  combined  effects  of  non-uniform  airflow  and  evaporator  superheat  were  examined 
by  blocking  the  upper  portion  of  the  test  evaporator  on  COIL-W  and  COIL-E  during 
cross-counter  flow  operation  for  wet  coil  conditions.  Figure  3. 4. 1.3.1  shows  an  idealized 
non-uniform  velocity  profile  for  a test  coil  with  the  upper  half  of  the  coil  partially 
blocked.  The  velocity  ratio  was  calculated  by  taking  the  average  of  the  1 5 velocity  points 
on  the  top  half  divided  by  the  average  of  the  1 5 velocity  points  on  the  lower  half  of  the 
test  evaporator.  Test  9 conditions  of  5.6  °C  (10.0  °F)  superheat  were  first  performed, 
then  the  blockage  was  applied  with  no  expansion  valve  adjustment  (test  9 A),  and  finally 
the  expansion  valves  were  adjusted  to  yield  5.6  °C  (10.0  °F)  superheat  on  all  circuits  (test 
9B).  During  these  tests  the  standard  airflow  rate  was  held  constant;  the  airflow  was  not 
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allowed  to  drop  when  the  blockage  was  added  regardless  of  the  significantly  higher 
pressure  drop.  Table  3.4. 1.3.1  shows  the  performance  of  the  coils  with  varying  degrees 
of  airflow  blockage. 


Figure  3 .4. 1 .3 . 1 : Idealized  velocity  profile  over  evaporator  with  upper  half  partially 
blocked 


Table  3.4.1 .3.1 : COIL-W  Performance  with  Non-Uniform  Airflow 


Velocity  Ratio 
(see  Figure  #) 

Test  name 

Test 

type' 

a 

o 

O 

Airflow, 
rrrVh  (scfm) 

Air-side 
Capacity,  W 
(Btu/h) 

Capacity 
Ratio, 
Q/Qlest  9 

Coil  Air 
Pressure 
Drop,  Pa 
(in  WG) 

1:1 

(Fig.  3.4.1. 3.2) 

W020522 

A 

9 

w 

1244  (732) 

6598  (22515) 

1 

37.9 

(0.152) 

1:1.5 

(Fig.  3.4.1. 3.3) 

W020523 

A 

9A 

w 

1252  (737) 

6351 (21670) 

0.96 

49.6 

(0.199) 

1:1.5 

W020524 

A 

9B 

w 

1249  (735) 

6535  (22298) 

0.99 

49.8 

(0.200) 

1:1 

W020528 

B 

9 

w 

1245  (733) 

6636  (22644) 

1 

38.1 

(0.153) 

1:2 

(Fig.  3.4. 1.3.4) 

W020528 

C 

9A 

w 

1239  (729) 

6179  (21085) 

0.93 

53.3 

(0.214) 

1:2 

W020529 

A 

9B 

w 

1237  (728) 

6307  (21521) 

0.95 

58.8 

(0.236) 

1)  Test  9:  uniform  airflow  with  superheat  on  all  circuits  of  5.6  °C  (10.0  °F),  9A:  same  expansion  valve  setting  as  test  9 but  with 
non-uniform  airflow  and  no  superheat  adjustment,  9B:  expansion  valves  adjusted  to  yield  5.6  °C  (10.0  °F)  superheat  on  all 
circuits  with  non-uniform  airflow. 
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Table  3.4. 1 .3.2:  COIL-E  Performance  with  Non-Uniform  Airflow 


Velocity  Ratio 
(see  Figure  #) 

Test  name 

Test 

type1 

Coil 

Airflow, 
m3/h  (scfm) 

Air-side 
Capacity,  W 
(Btu/h) 

Capacity 

Ratio, 

Q/Qlest  9 

Coil  Air 
Pressure 
Drop,  Pa 
(in  WG) 

1:1 

(Fig.  3.4.1.3.6) 

E020604 

A 

9 

E 

1276  (751) 

6985  (23833) 

1 

92.7 

(0.372) 

1:1.26 

(Fig.  3.4.1.3.7) 

E020604 

B 

9A 

E 

1281  (754) 

6933  (23655) 

0.99 

108.4 

(0.435) 

1:1.26 

E020605 

A 

9B 

E 

1281  (754) 

7029 (23984) 

1.01 

106.6 

(0.428) 

1:1 

E020607 

A 

9 

E 

1293  (761) 

6955  (23733) 

1 

84.7 

(0.340) 

1:1.36 

(Fig.  3.4. 1.3. 8) 

E020607 

B 

9A 

E 

1291  (760) 

6797  (23192) 

0.98 

102.9 

(0.413) 

1:1.36 

E020610 

A 

9B 

E 

1286  (757) 

6807  (23226) 

0.98 

101.9 

(0.409) 

1:1.62 

(Fig.  3.4.1.3.9) 

E020611 

A 

9A 

E 

1290(759) 

6751  (23034) 

0.97 

101.1 

(0.406) 

1:1.62 

E020612 

A 

9B 

E 

1274  (750) 

6914(23591) 

0.99 

101.4 

(0.407) 

1:1.75 

(Fig.  3.4.1.3.10) 

E020613 

A 

9A 

E 

1288  (758) 

6654  (22705) 

0.96 

105.4 

(0.423) 

1:1.75 

E020620 

A 

9B 

E 

1288  (758) 

6877  (23465) 

0.99 

103.4 

(0.415) 

1:2.59 

(Fig.  3.4.1.3.11) 

E020621 

A 

9A 

E 

1299  (764) 

6575  (22435) 

0.95 

98.9 

(0.397) 

1:2.59 

E020624 

A 

9B 

E 

1290  (759) 

6874  (23456) 

0.99 

101.9 

(0.409) 

] ) Test  9:  uniform  airflow  with  superheat  on  all  circuits  of  5.6  °C  (10.0  °F),  9A:  same  expansion  valve  setting  as  test  9 but  with 
non-uniform  airflow  and  no  superheat  adjustment,  9B:  expansion  valves  adjusted  to  yield  5.6  °C  (10.0  °F)  superheat  on  all 
circuits  with  non-uniform  airflow. 


Figure  3.4. 1.3. 2 shows  the  air  velocity  contour  map  for  COIL-W  with  no  obstructions 
present.  The  volumetric  flowrate  for  this  test  was  1244  m /h  (732  scfm)  with  an  average 
velocity  of  6437  m/h  (352  fpm)  and  standard  deviation  of  512  m/h  (28  fpm).  Any  non- 
uniformity in  the  unobstructed  evaporator’s  entrance  region  airflow  was  due  to  the 
dewpoint  sampling  tree,  thermocouple  grid,  and  fin  angles.  Figure  3.4. 1.3. 3 shows  the 
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non-uniform  velocity  contour  map  for  COIL-W  when  the  flow  was  obstructed  to  the 
upper  half  of  the  coil.  An  average  of  the  top  half  air  velocity  was  compared  to  the  bottom 
half  average  air  velocity  to  yield  the  velocity  ratio  of  1 to  1.5.  The  volumetric  flowrate 
for  this  test  was  1252  m3/h  (737  scfm)  with  a 4097  m/h  (224  fpm)  and  6163  m/h 
(337  fpm)  average  velocity  on  the  upper  and  lower  halves  of  the  coil,  respectively. 
Further  obstruction  was  added  to  produce  the  velocity  contours  seen  in  Figure  3. 4. 1.3. 4 at 
a velocity  ratio  of  1 to  2.  The  average  velocities  in  this  case  were  4005  m/h  (219  fpm) 
and  82 1 1 m/h  (449  fpm)  over  the  upper  and  lower  halves  of  the  coil,  respectively. 

The  imposed  airflow  blockage  in  the  case  of  the  1 to  1.5  velocity  ratio  would  have 
increase  fan  power  by  more  then  30  %.  For  the  1 to  2 velocity  ratio  case,  the  fan  power 
would  have  increased  by  at  least  54  % relative  to  the  uniform  airflow  case. 
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Figure  3. 4. 1.3. 2:  Uniform  airflow  velocity  (ft/min)  contour  map  for  COIL-W 


35 


Figure  3.4. 1.3. 3:  Non-uniform  velocity  (1:1.5)  (ft/min)  contour  map  for  COIL-W 
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Figure  3. 4. 1.3. 4:  Non-uniform  velocity  (1:2)  (ft/min)  contour  map  for  COIL-W 


Figure  3.4. 1.3. 5 shows  the  non-uniform  air  velocity  evaporator  capacity  relative  to  the 
uniform  air  velocity  capacity.  The  first  case  shows  the  effects  of  an  imposed  obstruction 
with  a fixed  area  expansion  device.  The  expansion  device  would  not  be  able  to  correct 
superheat  and  capacity  would  drop  by  4 % in  the  1:1.5  velocity  ratio  case  and  by  6.5  % in 
the  1 :2  velocity  ratio  case.  If  the  expansion  device  were  able  to  correct  each  circuit,  then 
the  performance  would  improve  closer  to  pre-obstruction,  but  for  both  cases  the  capacity 
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still  decreased  by  1.2  % and  4.5  %,  respectively.  The  COIL-W  non-uniform  velocity 
results  show  that  a more  non-uniform  velocity  ratio  produced  a higher  drop  in  capacity; 
and  superheat  correction  still  could  not  alleviate  the  entire  penalty. 


Figure  3.4.1. 3. 5:  Relative  capacities  for  COIL-W  non-uniform  airflow  (refer  to  Table 
3. 4. 1.3.1  for  a description  of  the  actual  tests) 

Figure  3.4. 1.3. 6 shows  the  air  velocity  contour  map  for  COIL-E  with  no  obstructions 
present.  The  volumetric  flowrate  for  this  test  was  1276  m/h  (751  scfm)  with  an  average 
velocity  of  6108  m/h  (334  fpm)  and  standard  deviation  of  512  m/h  (28  fpm).  Any  non- 
uniformity in  the  unobstructed  evaporator’s  entrance  region  airflow  was  due  to  the 
dewpoint  sampling  tree,  thermocouple  grid,  and  fin  angles.  Figure  3.4. 1.3. 7 shows  the 
non-uniform  velocity  contour  map  for  COIL-E  when  the  flow  was  obstructed  to  the  upper 
half  of  the  coil.  An  average  of  the  top  half  air  velocity  was  compared  to  the  bottom  half 
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average  air  velocity  to  yield  the  velocity  ratio  of  1 to  1 .26.  The  volumetric  flowrate  for 
this  test  was  1281  m3/h  (754  scfm)  with  a 5340  m/h  (292  fpm)  and  6748  m/h  (369  fpm) 
average  velocity  on  the  upper  and  lower  halves  of  the  coil,  respectively.  Further 
obstruction  was  added  to  produce  the  velocity  contours  seen  in  Figure  3.4. 1.3. 8 at  a 
velocity  ratio  of  1 to  1.36.  The  average  velocities  in  this  case  were  5340  m/h  (292  fpm) 
and  7279  m/h  (398  fpm)  over  the  upper  and  lower  halves  of  the  coil,  respectively.  More 
obstruction  was  added  to  produce  a 1 to  1.62  average  velocity  ratio  seen  in  Figure 
3.4. 1.3. 9.  Average  velocities  in  this  case  were  4609  m/h  (252  fpm)  and  7498  m/h 
(410  fpm)  over  the  upper  and  lower  halves  of  the  coil,  respectively.  Again,  obstruction 
was  added  to  produce  a 1 to  1.75  average  velocity  ratio  seen  in  Figure  3.4.1.3.10. 
Average  velocities  in  this  case  were  3621  m/h  (198  fpm)  and  6346  m/h  (347  fpm)  over 
the  upper  and  lower  halves  of  the  coil,  respectively.  The  final  obstruction  was  then  added 
to  produce  a velocity  ratio  of  1 to  2.59.  Average  velocities  in  this  case  were  3219  m/h 
(176  fpm)  and  8358  m/h  (457  fpm)  over  the  upper  and  lower  halves  of  the  coil, 


respectively. 
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Figure  3.4. 1.3. 7:  Non-uniform  velocity  (1:1.26)  (ft/min)  contour  map  for  COIL-E 


41 


Figure  3.4.1. 3. 8:  Non-uniform  velocity  (1:1.36)  (ft/min)  contour  map  for  COIL-E 
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Figure  3.4. 1.3. 9:  Non-uniform  velocity  (1:1.62)  (ft/min)  contour  map  for  COIL-E 
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Figure  3.4.1.3.10:  Non-uniform  velocity  (1:1.75)  (ft/min)  contour  map  for  COIL-E 
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Figure  3. 4. 1.3.1 1:  Non-uniform  velocity  (1:2.59)  (ft/min)  contour  map  for  COIL-E 


Figure  3.4.1.3.12  shows  the  non-uniform  air  velocity  evaporator  capacity  relative  to  the 
uniform  air  velocity  capacity  for  all  COIL-E  tests  represented  above  and  shown  in  Table 
3.4. 1.3.2.  The  cases  labeled  “Test  9A”  show  the  effects  of  an  imposed  obstruction  with  a 
fixed  area  expansion  device.  Cases  labeled  “Test  9B”  have  the  superheat  adjusted  back 
to  5.6  °C  (10.0  °F)  on  each  circuit.  The  first  case  had  a velocity  ratio  of  1 to  1.26  and  a 
capacity  specific  airflow  of  183  m3/kWh  (378  scffn/ton).  Less  than  1 % in  capacity  was 
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lost  due  to  the  obstruction,  and  this  lost  capacity  was  recovered  when  superheat  was 
corrected  (capacity  was  corrected  within  the  uncertainty  of  our  measurement  at 
approximately  3.5  %).  At  a higher  absolute  airflow  of  186  m3/kW  h (385  scfm/ton),  the 
losses  in  capacity  with  the  obstruction  were  greater. 

The  losses  in  capacity  with  no  superheat  correction  ranged  from  slightly  more  than  2 % at 
the  1 to  1.36  velocity  ratio  to  approximately  5.5  % at  the  1 to  2.59  velocity  ratio.  When 
the  superheat  was  corrected,  much  of  the  loss  in  capacity  was  recovered.  Please  note  that 
all  of  these  tests  were  performed  at  a constant  airflow  rate.  As  shown  in  Table  3.4. 1.3.1 
and  3. 4. 1.3. 2,  the  air  pressure  drop  across  the  evaporator  increased  substantially  when  the 
blockage  was  imposed.  This  would  translate  into  much  higher  fan  power  requirements 
and  a subsequently  lower  COP.  In  the  case  of  the  1 to  1 .36  velocity  profile,  fan  power 
would  have  increased  by  21  % with  the  imposed  blockage  (power  equals  flowrate  times 
the  pressure  drop).  For  the  1:1.62  velocity  profile,  fan  power  would  have  increased  by 
19  %.  For  the  1:1.75  velocity  ratio,  fan  power  would  have  increased  by  24  %.  For  the 
highest  blockage  test  with  a velocity  ratio  of  1:2.59,  the  fan  power  would  have  increased 
by  at  least  20  %. 
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Capacity  Relative  to  Uniform  Airflow 


183  m /kWh  (378  scfm/ton) 


Test9A  Test9B  Test9A  Test9B  Test9A  Test9B  Test9A  Test9B  Test9A  Test9B 


Figure  3.4.1.3.12:  Capacity  relative  to  uniform  airflow  for  COIL-E  (refer  to  Table 

3. 4. 1.3. 2 for  a description  of  the  tests) 
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3.4.2  Cross-Parallel  Air/Refrigerant  Flow  Configuration  Tests 


Table  3.4.2. 1 shows  the  performance  of  COIL-W  and  COIL-E  with  refrigerant  circuited 
in  cross-counter  flow  and  cross-parallel  flow  with  a capacity  specific  airflow  rate  of 
1 93  m3/h  (400  scfm/ton).  COIL-W  airflow  rates  for  parallel  flow  and  counter  flow  were 
912  m3/h  (537  scfm)  and  1240  m3/h  (730  scfm),  respectively.  COIL-E  airflow  rates  for 
parallel  flow  and  counter  flow  were  895  m3/h  (527  scfm)  and  1288  m3/h  (758  scfm), 
respectively.  The  main  information  to  be  gained  from  this  comparison  is  the  rapid 
reduction  in  capacity  that  follows  any  increase  in  superheat  during  parallel  flow  (Figure 
3.4.2. 1).  This  was  evident  for  both  coils  tested  under  parallel  flow  conditions. 

As  the  superheat  increased  from  5.6  °C  (10.0  °F)  to  16.7  °C  (30.0  °F),  the  refrigerant  and 
air  temperatures  became  pinched  very  quickly.  During  all  tests  the  evaporator  pressure 
was  fixed  to  give  an  evaporating  saturation  temperature  of  7.2  °C  (45.0  °F).  When 

superheat  was  increased  to  16.7  °C  (30.0  °F),  the  exiting  refrigerant  temperature 
approached  23.9  °C  (75.0  °F);  pinching  of  the  two  streams  occurred  rapidly.  This  is 
evident  from  examining  the  capacity  of  test  10  in  parallel  flow. 

Figure  3.4.2. 1 shows  that  COIL-W  capacity  in  parallel  flow  decreased  by  84.8  % as  the 
superheat  increased  from  5.6  °C  (10.0  °F)  to  16.7  °C  (30.0  °F);  test  9 to  test  10.  Tests  1 1 
and  12  also  produced  very  low  capacity  due  to  two  of  the  three  circuits  having  a 16.7  °C 
(30.0  °F)  superheat.  Although  tests  with  the  center  circuit  flooding  produced  higher 
capacity  than  the  top  circuit  flooding  for  both  COIL-W  and  COIL-E. 
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Table  3.4.2. 1:  Counter  and  Parallel  Flow  Performance  of  COIL-W  and  COIL-E 


Test  Name 

Test  # 

Coil  Designation 

Volumetric  Flowrate  of 
Air  m3/h  (cfm) 

Coil 

Surface 

Overall  Superheat 

Superheats  in  Individual  Circuits 

145-Q,a 

(300-Q) 

193-Qla 

(400-Q) 

242-Qla 

(500-Q) 

d 

-4—> 

o 

£ 

5.6  °C 
(10.0  °F) 

16.7  °C 
(30.0  °F) 

5.6  °C 
(10.0  °F) 

5.6  °C(10.0  °F) 

5. 6/5. 6/5. 6 
(10/10/10) 

16.7/16.7/16.7 

(30/30/30) 

16.7/*/16.7 

(30/*/30) 

’•‘/l  6.7/1 6.7 
(*/30/30) 

Qtest 

Qtest/ 

Qlb 

Qtest 

Qtest/ 

Qlb 

Qtest  r 

)tes(/ 1 

Qtest 

5test/ 

3lb 

W (Btu/h) 

W (Btu/h) 

W (Btu/h)  Q 

yb  \ 

V (Btu/h)  ( 

Cross-Counter  Air/Refrigerant  Flow 

W020207B 

_ 

9 

W 

X 

X 

6507 

(22203) 

l 

I 

W020530A 

10 

w 

X 

X 

3722 

(12701) 

0.57 

W020531 A 

11 

W 

X 

X 

3837 

(13091) 

0.59 

W020215B 

12 

W 

X 

X 

3830 

(13067) 

0.59 

E020607A 

9 

E 

X 

X 

6955 

(23733) 

l 

W020318A 

10 

E 

X 

X 

4865 

(16599) 

0.70 

W020318B 

1 1 

E 

X 

X 

5485 

(18715) 

0.79 

W020319A 

12 

E 

X 

X 

4735 

(16157) 

0.68 

Cross-Parallel  Air/Refrigerant  Flow 

W020304A 

9 

W 

X 

X 

4732 

(16146) 

, 

W020311B 

10 

W 

X 

X 

721 

(2461) 

0.15 

W020306A 

1 1 

W 

X 

X 

2605 

(8887) 

0.55 

W020307A 

12 

W 

X 

X 

2143 

(7311) 

0.45 

E020403A 

9 

E 

X 

X 

4549 

(15523) 

l 

E020404A 

10 

E 

X 

X 

1017 

(3470) 

0.22 

E020408A 

11 

E 

X 

X 

3373 

(11508) 

0.74 

E020409A 

12 

E 

X 

X 

2797 

(9543) 

0.61 

1 a : capacity  determined  at  1 93  nr’/h  (400  scfm/ton) 
lb  : capacity  at  test  9 for  the  coil  specified. 
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Figure  3.4.2. 1:  Cross-parallel  flow  capacity  comparison  for  COIL-W  and  COIL-E 
relative  to  their  performance  at  test  9 with  cross-counter  flow 

COIL-E  capacity  in  parallel  flow  was  3.9  % less  than  COIL-W  at  the  test  9 condition  of 
193  m3/kWh  (400  scfm/ton).  This  is  the  opposite  of  the  counter  flow  capacity  results 
where  COIL-E  had  a 5.4  % greater  capacity  than  COIL-W.  COIL-E  produced  a slightly 
lower  capacity  than  COIL-W  at  the  test  9 parallel  flow  conditions.  This  was  mainly  due 
to  the  196  m3/kWh  (406  scfm/ton)  for  COIL-W  and  the  193  nvVkWh  (399  scfm/ton)  for 
COIL-E.  The  ideal  airflow  would  have  produced  193  m3/kWh  (400  scfm/ton)  for  both 
coils,  but  COIL-W  airflow  was  high  while  COIL-E  airflow  was  slightly  low.  These 
differences  in  airflow  produced  the  accompanying  difference  in  capacity.  Also,  a 
secondary  factor  in  the  lower  capacity  of  COIL-E  than  COIL-W  may  be  due  to  the  more 
rapid  pinching  of  COIL-E  than  COIL-W  due  to  the  higher  air-side  heat  transfer  of  the 
enhanced  fins. 
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COIL-E  capacity  in  parallel  flow  decreased  by  77.6  % as  the  superheat  increased  from 
5.6  °C  (10.0  °F)  to  16.7  °C  (30.0  °F);  test  9 to  test  10.  Again,  rapid  pinching  of  the  two 
fluid  streams  reduced  capacity. 
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4 EVAPORATOR  MODELING  AND  SIMULATIONS 

4.1  Background  on  Evaporator  Model  EVAP5 

Modeling  of  finned  tube  heat  exchangers  started  at  NIST  with  a tube-by-tube  simulation 
model  originally  formulated  by  Chi  (1979).  Over  the  years,  the  model  underwent 
significant  upgrades,  which  are  documented  in  Domanski  and  Didion  (1983),  Domanski 
(1991),  Lee  and  Domanski  (1997),  and  Domanski  (1999b).  These  upgrades  included  the 
capability  to  account  for  air  maldistribution  and  its  interaction  with  refrigerant 
distribution,  the  extension  to  zeotropic  mixtures,  the  extension  to  new  refrigerant 
property  representations,  and  implementations  of  new  simulation  correlations.  The  1 999 
upgrade  equipped  the  evaporator  model  with  a graphical  user  interface  (GUI).  The  GUI 
serves  as  a pre-  and  post-processor;  it  facilitates  preparation  of  simulation  input  data, 
including  the  layout  of  refrigerant  circuitry,  and  allows  the  user  to  display  detailed 
performance  results  for  individual  tubes  after  a simulation  run  is  completed.  In  2002 
under  a parallel  ARTI-21CR/605-500100  project  (Domanski  and  Payne,  2002),  the 
EVAP-COND  simulation  package  was  developed  that  included  the  evaporator  and 
condenser  models,  EVAP5  and  COND5,  working  under  the  same  GUI.  Both  models 
were  validated  with  experimental  data  taken  with  R22  and  R410A  heat  exchangers  at 
NIST.  Since  the  current  project  and  ARTI-21CR/605-500100  project  partially 
overlapped  in  time,  the  EVAP-COND  attached  to  the  ARTI-21CR/605-500100  report 
included  some  of  upgrades  developed  under  this  study,  e.g.  the  option  to  simulate  the 
evaporator  together  with  a refrigerant  distributor.  Not  included  in  the  release  version  of 
EVAP-COND  are  the  upgrades  to  the  evaporator  model  that  account  for  longitudinal  fin 
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conduction,  which  were  needed  to  perform  simulations  with  controlled  refrigerant 
superheats  at  individual  refrigerant  circuits. 

4.2  Description  of  EVAP5 

4.2.1  Modeling  Approach 

Figure  4. 2. 1.1  presents  the  refrigerant  circuitry  and  air  velocity  representation  used  by 
EVAP5.  The  tube-by-tube  modeling  approach  recognizes  each  tube  as  a separate  entity 
for  which  the  model  performs  simulation  calculations.  These  calculations  are  based  on 
inlet  refrigerant  and  air  parameters,  their  properties,  and  mass  flow  rates.  The  simulation 
begins  with  the  inlet  refrigerant  tubes  and  proceeds  to  successive  tubes  along  the 
refrigerant  path  in  the  heat  exchanger.  At  the  outset  of  the  simulation,  the  air  temperature 
is  only  known  for  the  tubes  in  the  first  row  and  has  to  be  estimated  for  the  remaining 
tubes.  A successful  simulation  run  requires  several  passes  (iterations)  through  the 
refrigerant  circuitry,  each  time  updating  inlet  air  and  refrigerant  parameters  for  each  tube. 
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Fins 


Air  velocity  profile 


Coil  height 


Figure  4.2. 1 . 1 : Representation  of  air  distribution  and  refrigerant  circuitry  in  EVAP5 


Heat  transfer  calculations  start  by  calculating  the  heat-transfer  effectiveness,  £,  by  one  of 
the  applicable  relations  (Kays  and  London,  1984).  With  the  air  temperature  changing  due 
to  heat  transfer,  the  selection  of  the  appropriate  relation  for  £ depends  on  whether  the 
refrigerant  undergoes  a temperature  change  during  heat  transfer.  Once  £ is  determined, 
heat  transfer  from  air  to  refrigerant  is  obtained  using  equation  4.2. 1 . 1 . 


Qa  =maCpa(Tai  -Tn)e 


(4.2. 1.1) 


The  overall  heat-transfer  coefficient,  U,  is  calculated  by  equation  4.2. 1.2,  which  sums  up 
the  individual  heat-transfer  resistances  between  the  refrigerant  and  the  air. 
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(4.2. 1.2) 


The  first  term  of  equation  4.2. 1.2  represents  the  refrigerant-side  convective  resistance. 
The  second  term  is  the  conduction  heat-transfer  resistance  through  the  tube  wall,  and  the 
third  term  accounts  for  the  conduction  resistance  through  the  water  layer  on  the  fin  and 
tube.  The  fourth  term  represents  the  contact  resistance  between  the  outside  tube  surface 
and  the  fin  collar.  The  fifth  term  is  the  convective  resistance  on  the  air-side  where  the 
multiplier  (1+a)  in  the  denominator  accounts  for  the  latent  heat  transfer  on  the  outside 
surface.  For  a dry  tube  a = 0.0  and  1 /hi  = 0.0.  Once  the  heat  transfer  rate  from  the  air  to 
the  refrigerant  is  calculated,  the  tube  wall  and  fin  surface  temperatures  can  be  calculated 
directly  using  heat-transfer  resistances.  Then,  the  humidity  ratios  for  the  saturated  air  at 
the  wall  and  fin  temperatures  are  calculated,  and  mass  transfer  from  the  moist  air  to  the 
tube  and  fin  surfaces  is  determined.  For  more  detailed  information  on  heat  transfer 
calculations  refer  to  Domanski  (1991). 


4.2.2  Heat  Transfer  and  Pressure  Drop  Correlations 

EVAP5  uses  the  following  correlations  for  calculating  heat  transfer  and  pressure  drop. 

Air  Side 

• heat-transfer  coefficient  for  flat  fins:  Wang  et  al.  (2000) 

• heat-transfer  coefficient  for  wavy  fins:  Wang  et  al.  (1999a) 

• heat-transfer  coefficient  for  slit  fins:  Wang  et  al.  (2001) 
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• heat-transfer  coefficient  for  louver  fins:  Wang  et  al.  (1999b) 

• fin  efficiency:  Schmidt  method,  described  in  McQuiston  et  al.  (1982) 

A correlation  for  calculating  the  tube-collar  junction  resistance  is  not  listed  because  all 
air-side  heat  transfer  correlations  authored  by  Wang  include  the  heat  transfer  resistance 
between  the  tube  and  collar. 


Refrigerant  Side 

• single-phase  heat-transfer  coefficient,  smooth  tube:  McAdams,  described  in 
ASHRAE  (2001) 

• evaporation  heat-transfer  coefficient  up  to  80%  quality,  smooth  tube:  Jung  and 
Didion  (1989) 

• evaporation  heat-transfer  coefficient  up  to  80%  quality,  rifled  tube:  Jung  and 
Didion  (1989)  correlation  with  a 1 .9  enhancement  multiplier  suggested  by 
Schlager  et  al.  (1989) 

• mist  flow,  smooth  and  rifled  tubes:  linear  interpolation  between  heat  transfer 
coefficient  values  for  80  % and  1 00  % quality 

• single-phase  pressure  drop,  smooth  tube:  Petukhov  (1970) 

• two-phase  pressure  drop,  smooth  tube,  lubricant-free  refrigerant:  Pierre  (1964) 

• two-phase  pressure  drop,  rifled  tube:  Pierre  (1964)  correlation  for  smooth  tube 
with  a 1.4  multiplier  suggested  by  Schlager  et  al.  (1989) 

• single-phase  pressure  drop,  return  bend,  smooth  tube:  White,  described  in 
Schlichting  (1968) 

• two-phase  pressure  drop,  return  bend,  smooth  tube:  Chisholm,  described  in 
Bergles  et  al.  (1981).  The  length  of  a return  bend  depends  on  the  relative 
locations  of  the  tubes  connected  by  the  bend.  This  length  is  accounted  for  in 
pressure  drop  calculations. 
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4.2.3  Representation  of  Refrigerant  Properties 


Representation  of  thermodynamic  and  transport  properties  is  based  on  REFPROP6 
property  routines  (McLinden  et  al.,  1998).  Because  EVAP5  simulations  are 
computationally  intensive,  using  a refrigerant  property  look-up  tables  is  a practical 
necessity  if  simulation  runs  are  expected  to  take  less  than  60  seconds.  This  is  particularly 
true  in  case  of  REFPROP6  for  which  property  calculations  are  several  times  more  CPU 
demanding  than  for  REFPROP5.  EVAP-COND  employs  a pressure-enthalpy-based 
system  of  look-up  tables,  which  includes  eight  different  routines  that  retrieve  the  desired 
state  or  transport  property.  The  look-up  scheme  is  applicable  to  single  component 
refrigerants  and  refrigerant  mixtures.  If  a given  refrigerant  state  falls  outside  the  range  of 
the  look-up  table,  then  EVAP-COND  calls  a REFPROP6  refrigerant  property  routine 
directly. 

Since  REFPROP6  property  calculations  may  not  converge  occasionally  , particularly 
during  phase  equilibrium  calculations  for  refrigerant  mixtures,  EVAP-COND  employs  an 
error  evasive  scheme.  Under  this  scheme,  EVAP-COND  attempts  to  obtain  a given 
property  even  if  REFPROP  flash  calculations  do  not  converge,  e.g.,  if  a routine  PHFLSH 
crashes,  a routine  that  uses  TPRHO  is  invoked  to  attempt  to  iteratively  match  TPRHO’s 
enthalpy  value  with  the  known  (target)  value.  If  both  REFPROP  flash  calculations  do  not 
converge,  then  the  data  point  in  the  refrigerant  look-up  table  is  flagged  and  look-up  table 
routines  iterate  properties  for  this  point  using  refrigerant  properties  in  the  neighboring 
nodes  of  the  table. 
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4.3  Modeling  Issues 

The  following  four  sections  present  four  major  modeling  issues  that  received  special 
consideration  during  this  study. 

4.3.1  Refrigerant  Distribution 

Simulating  refrigerant  distribution  is  an  important  aspect  of  heat  exchanger  simulation 
because  of  its  impact  on  the  heat  exchanger  performance.  It  is  also  know  that  in  some 
designs  a non-uniform  air  distribution  may  affect  refrigerant  distribution.  In  a heat 
exchanger  with  multiple  circuits,  refrigerant  distributes  itself  in  appropriate  proportions 
so  that  the  refrigerant  pressure  drop  from  the  inlet  to  the  outlet  for  all  circuits  is  the  same. 
In  the  context  of  simulating  refrigerant  distribution,  a refrigerant  circuit  starts  at  the  point 
of  the  first  split  of  refrigerant  stream  after  leaving  the  condenser  and  ends  at  the  final 
merging  point  before  entering  the  suction  line  leading  the  refrigerant  to  the  compressor. 
If  the  refrigerant  enters  the  evaporator  by  a single  tube,  the  first  split,  if  any,  will  exist 
within  the  coil  assembly.  If  the  evaporator  has  several  inlet  tubes  and  a refrigerant 
distributor  is  used,  the  first  refrigerant  split  typically  occurs  at  the  inlet  to  the  distributor 
tubes  just  after  the  expansion  process  in  a thermostatic  expansion  valve  (TXV)  or  a short 
tube  restrictor.  Note,  that  in  this  design,  refrigerant  pressures  and  temperatures  at 
different  inlet  tubes  may  be  different,  as  graphically  shown  in  Figure  4. 3. 1.1.  Such 
different  refrigerant  pressure  and  temperature  profiles  also  occurred  during  the  tests  with 
controlled  uneven  exit  superheats  (refrigerant  distributions),  namely  tests  3,  4,  7,  8,  11, 
and  12. 
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Figure  4.3. 1 . 1 : Possible  refrigerant  pressure  profiles  in  a three-circuit  evaporator  fed  by  a 
refrigerant  distributor 

Under  this  project,  two  simulation  methods  were  developed  to  simulate  evaporator 
performance  with  controlled  refrigerant  superheats  at  the  evaporator  outlet  tubes.  The 
first  scheme  involving  a general  model  for  a refrigerant  distributor  was  introduced  into 
EVAP-COND  as  one  of  the  eight  evaporator  simulation  options.  When  the  evaporator  is 
simulated  using  this  option,  the  refrigerant  operating  input  data  are  the  condenser  exit 
bubble-point  temperature,  condenser  subcooling,  evaporator  exit  dew-point  temperature, 
and  evaporator  superheat  (in  the  release  version  of  EVAP-COND,  condenser  subcooling 
and  evaporator  superheat  are  imposed  as  5.0  °C  (9.0  °F)).  For  this  option,  as  the  first  task 
EVAP-COND  runs  preliminary  simulations  to  establish  dimensions  of  the  refrigerant 
distributor  tubes  that  would  inflict  a 70  kPa  (10.2  lb/in2)  refrigerant  pressure  drop.  Once 
the  distributor  tubes  are  sized,  EVAP-COND  proceeds  to  main  simulations  in  which  it 


59 


establishes  refrigerant  distribution  between  different  circuits  based  on  the  total  pressure 
drop.  This  total  pressure  drop  includes  the  pressure  drop  in  a given  distributor  tube  and 
the  refrigerant  circuit  in  the  coil  assembly  it  feeds.  In  the  test  version  used  in  this  study, 
EVAP-COND  additionally  solicits  a “restriction  factor”  for  each  distributor  tube,  which 
acts  as  a multiplier  to  the  pressure  drop  calculated  by  the  program.  By  inputting  values 
different  from  1 .0,  the  user  can  control  refrigerant  distribution  and  refrigerant  superheat 
at  different  evaporator  exit  tubes.  The  program  iterates  the  refrigerant  mass  flow  rate 
until  the  overall  superheat  is  reached  at  the  evaporator  exit.  Figures  4. 3. 1.2  and  4. 3. 1.3 
present  the  eight  refrigerant  input  data  options  and  the  input  data  window  for  EVAP- 
COND  simulations  involving  a refrigerant  distributor,  respectively. 


Outlet  pressure  and  siperheat 
Outlet  pressure  and  quality 

Outlet  sat  temp,  and  superheat 
Outlet  sat.  temp,  and  quality 

Inlet  pressure  and  quality 
Inlet  sat.  temp,  and  quality 


v*  Cond.  bubble  pt:  and  evap.  sat.  temp. 


Outlet  sat.  temp,  and  min.  superheat 


Figure  4.3.1 .2:  EVAP-COND  refrigerant  input  data  options  for  evaporator  simulations 


Figure  4.3.1 .3:  EVAP-COND  input  data  window  for  simulations  involving  a refrigerant 
distributor 
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While  the  simulation  option  involving  a refrigerant  distributor  is  useful  for  a coil  designer 
for  typical  simulations,  this  option  proved  to  be  somewhat  impractical  for  controlling 
individual  exit  superheats  when  we  tried  to  reproduce  our  test  results.  This  was  due  to  a 
non-linear  dependence  of  individual  exit  superheats  on  the  “restriction  factors”.  For  this 
reasonanother  simulation  scheme  was  developed  in  which  the  user  directly  assigns 
refrigerant  distribution  between  different  circuits.  The  operating  conditions  are  as  shown 
in  Figure  4. 3. 1.4  with  the  refrigerant  inlet  quality  and  distribution  (in  fractions)  solicited 
by  the  follow-up  DOS  window.  While  holding  the  refrigerant  distribution  constant,  the 
program  iterates  the  overall  refrigerant  mass  flow  rate  and  inlet  pressures  at  individual 
inlet  tubes  to  converge  on  the  target  exit  pressure  (the  same  for  each  exit  tube)  and 
overall  target  superheat.  Different  individual  superheats  can  be  obtained  by  specifying 
different  refrigerant  distributions.  Eventually,  all  simulation  results  for  this  study  were 
obtained  using  the  second  scheme. 


Figure  4. 3. 1.4:  EVAP-COND  input  data  window  for  simulations  with  specified  overall 
evaporator  exit  saturation  temperature  and  superheat. 
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4.3.2.  Air-side  Heat  Transfer  Correlations 


Often  the  most  significant  part  of  heat  transfer  resistance  between  the  air  and  refrigerant 
is  on  the  air-side  of  the  heat  exchanger.  For  this  reason,  a literature  review  on  the  latest 
air-side  heat  transfer  correlations  was  performed  at  the  outset  of  this  study.  Of  our 
particular  interest  were  correlations  for  wavy  and  lanced  (slit)  fins  - the  two  fin  types 
used  in  the  evaporators  tested  under  this  project. 

Figure  4.3.2. 1 compares  the  predictions  of  different  correlations  available  in  the 
literature.  These  predictions  were  calculated  for  typical  fin  designs  for  a three-depth-row 
heat  exchanger.  The  layout  of  different  prediction  lines  in  the  figure  may  serve  as  an 
explanation  why  predicting  performance  of  a finned-tube  heat  exchanger  may  be 
difficult.  For  wavy  fins,  the  correlation  by  Wang  et  al.  (1999a)  and  Kim  et  al.  (1997)  are 
in  a close  agreement,  while  the  correlation  by  Webb  (1990)  calculates  heat  transfer 
coefficients  up  to  50  % lower  that  the  two  first  methods.  In  the  air  velocity  range  of  1.8 
m/s  (5.9  ft/s)  to  2.1  m/s  (6.9  ft/s),  the  Webb  correlation  breaks  sharply  due  to  switching 
between  two  different  algorithms  with  a changing  air-side  Reynolds  number.  At  some 
point  the  Webb  correlations  provide  a value  for  the  wavy  fin  heat  transfer  coefficient  that 
is  lower  than  that  for  a flat  fin,  which  is  not  a realistic  prediction. 

For  slit  (lanced)  fins,  the  correlations  by  Nakamura  and  Xu  (1983)  and  Wang  et  al. 
(2001)  may  differ  by  more  than  40  %,  depending  on  air  velocity.  This  spread  may  be 
indicative  of  the  general  fact  that  some  correlations  do  not  predict  well  outside  the 
geometries  for  which  they  were  developed.  A measurement  uncertainty  in  one  or  both 
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experiments  may  also  be  a contributing  factor  to  this  large  discrepancy.  In  addition,  it 
should  be  noted  that  the  Nakayama  and  Xu  (1983)  predictions  do  not  approach  zero  at  air 
velocities  below  2 m/s  (6.5  ft/s),  the  trend  exhibited  by  the  other  correlations.  Regarding 
louver  fins,  the  correlation  by  Wang  et  al.  (1999b)  shows  a step  change  in  the  1.5  m/s 
(4.9  ft/s)  to  1.8  m/s  (5.9  ft/s)  range  caused  by  using  two  different  algorithms,  similar  to 
the  Webb  correlation  for  the  wavy  fins. 

The  analysis  of  relative  predictions  of  the  air-side  heat  transfer  coefficient  provided  the 
reason  for  replacing  the  existing  correlations  in  EVAP5  with  correlations  published  by 
Wang  and  his  co-workers  for  all  types  of  fins,  i.e.,  flat,  wavy,  louver,  and  slit  fins.  It  was 
judged  that  a better  degree  of  prediction  consistency  can  be  obtained  with  all  correlations 
developed  by  the  same  author.  Still,  the  reader  should  note  a reservation  regarding  the 
louver  fin  correlation,  which  did  not  provide  smooth  predictions  in  Figure  4.3.2. 1. 

In  conclusion,  we  have  to  recognize  the  spread  in  performance  between  different 
enhanced  fins,  either  realistic  or  perhaps,  in  some  instances,  overstated  by  correlations. 
To  accommodate  these  differences  and  facilitate  accurate  evaporator  model  predictions, 
EVAP-COND  provides  an  option  that  allows  the  user  to  “tune”  evaporator  simulated 
performance  to  the  laboratory  data  by  specifying  a “correcting  parameter”  for  the  air-side 
heat  transfer  coefficient  (such  correcting  parameters  are  also  allowed  for  the  refrigerant- 
side  heat  transfer  coefficient,  and  refrigerant  pressure  drop). 
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Heat  Transfer  Coefficient  (W/m2K) 


2 3 4 

Velocity  (m/s) 


Figure  4.3.2. 1 : Comparison  of  air-side  heat  transfer  correlations 
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4.3.3.  Internal  Heat  Transfer  in  a Finned-Tube  Evaporator 


The  current  study  stipulated  evaporator  tests  with  a superheat  as  high  as  16.7  °C 
(30.0  °F)  to  assess  capacity  degradation  at  large  and  uneven  superheats.  As  presented  in 
previous  sections,  these  tests  for  COIL-E  and  COIL-EC  produced  rather  interesting  data 
suggesting  that  internal  heat  transfer  within  the  evaporator  metal  body  may  be  the  culprit 
for  significant  capacity  degradation.  Figure  4.3.3. 1 presents  measured  capacities  at 
193  m3/kWh  (400  cfm/ton)  and  includes  capacities  predicted  by  EVAP5  (internal  heat 
transfer  not  considered).  The  figure  shows  that  tested  capacities  at  5.6  °C  (10  °F)  overlap 
for  both  coils  and  are  reasonably  well  predicted  by  EVAP5.  However,  at  16.7  °C  (30  °F) 
superheat  COIL-E  capacity  was  tested  to  be  346  W (1180  Btu/h)  less  that  COIL-EC 
capacity  and  966  W (3295  Btu/h)  less  than  the  simulated  value.  The  lower  capacity 
degradation  for  COIL-EC  can  be  explained  by  the  fin  cuts,  which  physically  separated 
different  tube  depth  rows  and  disallowed  heat  transfer  between  neighboring  rows. 
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Figure  4.3.3. 1:  Tested  and  predicted  capacities  for  COIL-E  and  COIL-EC  at  5.6  °C  (10 
°F)  and  16.7  °C  (30  °F)  superheats  (Test  9 and  Test  10  operating 
conditions.  Internal  heat  transfer  within  the  evaporator  metal  body  not 
considered.) 

The  following  two  sections  discuss  the  two  modes  of  the  internal  heat  transfer, 
longitudinal  tube  conduction  and  longitudinal  fin  conduction  (tube-to-tube  heat  transfer), 
and  their  impact  on  evaporator  performance. 


4.3.3. 1 Longitudinal  Tube  Conduction 

The  general  theory  states  that  if  a temperature  gradient  exists  in  a wall  of  a heat 
exchanger,  then  conduction  heat  transfer  will  occur  along  that  wall  and  may  therefore 
degrade  the  performance  of  the  heat  exchanger.  Kays  and  London  (1984)  identified  the 
major  parameters  affecting  the  magnitude  of  the  performance  degradation  due  to  this 
phenomenon  as  follows: 
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The  magnitude  of  the  performance  degradation  becomes  larger  with  increasing  X, 
Cmin/Cmax,  and  NTU.  Kays  and  London  stated  that  this  reduction  in  performance  is  seen 
in  heat  exchangers  designed  for  high  effectiveness  (£>0.9),  however,  they  did  not  provide 
much  of  a quantitative  analysis.  Ranganayakulu  et.  al.  (1996)  carried  out  a series  of  finite 
element  simulations  to  quantify  the  magnitude  of  the  performance  degradation  in  a heat 
exchanger  due  to  longitudinal  heat  conduction..  The  results  of  their  simulations  are 
represented  by  the  “conduction  effect  factor”,  x,  in  terms  of  the  effectiveness  with  no 
longitudinal  conduction  effects,  £NC,  and  the  effectiveness  considering  longitudinal 

conduction  effects,  £wc . 
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NC 


(4.3.3. 1.2) 


The  conduction  effect  factor  can  be  read  from  the  charts  presented  in  the  paper  for  given 
£,A,  Cmin/Cmax,  and  NTU.  Ranganayakulu  et.  al.  suggested  0.8  as  the  effectiveness  limit 
below  which  the  impact  of  longitudinal  conduction  is  negligible. 


With  the  theory  and  the  results  from  the  numerical  simulations  at  hand,  the  impact  of 
longitudinal  tube  conduction  for  a typical  finned-tube  evaporator  was  examined.  Using 
EVAP5,  a 10.6  kW  (3  ton)  evaporator  was  simulated  to  identify  the  tubes  with  two-phase 
R-22  (in  which  the  longitudinal  heat  conduction  does  not  occur)  and  the  tubes  with  a 
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superheated  refrigerant  (in  which  longitudinal  heat  conduction  does  take  place).  Then  the 
capacity  penalty  for  the  superheated  tubes  was  calculated  as  a fraction  of  capacity  of 
these  tubes  and  as  a fraction  of  the  evaporator  capacity. 


Figure  4. 3. 3. 1.1  displays  an  evaporator  side  view  with  a schematic  of  the  refrigerant 
circuitry.  Figure  4. .3. 3. 1.2  contains  a Coil  Design  Data  window  from  EVAP-COND  with 
the  coil  design  information,  and  Figure  4.3. 3. 1.3  presents  the  operating  conditions  of  the 
evaporator. 


Figure  4. 3. 3.1 .1 : Refrigerant  circuitry  configuration  for  the  analyzed  R-22  evaporator 
(inlet  tube:  tube  # 24,  outlet  tubes:  tube  # 1 and  tube  # 16;  ) 
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rData  for  a section  - 
No.  of  tubes  in  depth  row  #1 : JTb 
No.  of  tubes  in  depth  row  8 2. 

No.  of  tubes  in  depth  row  83: 

No.  of  tubes  in  depth  row  84: 

No.  of  tubes  in  depth  row  85: 


rTtie  data- 


. 
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0.386 
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Thickness  1 mm 
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Wavy  ?r 

Thermal  conductivity  | kW/(m.C) 

0.2216 

■ 
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> > 
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Figure  4. 3. 3. 1.2:  EVAP-COND  window  with  evaporator  design  information 


Figure  4.3.3. 1 .3:  EVAP-COND  window  with  evaporator  operating  conditions 


For  the  evaporator  exit  superheat  of  8.0  °C  (14.4  °F),  the  simulations  showed  that  only  5 
of  the  48  tubes  in  the  heat  exchanger  have  superheated  refrigerant  and  experience  a 
temperature  change.  This  means  that  43  of  the  tube  passes  in  the  heat  exchanger  will  not 
experience  any  axial  heat  conduction  because  there  will  be  no  temperature  difference  for 
this  to  occur  (neglecting  the  marginal  drop  of  saturation  temperature  due  to  the  pressure 
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drop).  An  example  of  a tube  with  superheated  vapor,  tube  number  15,  has  the  following 
values  for  the  aforementioned  parameters. 


( w ) 
386 

(l,1932E-5m2) 

_ kAw  _ 

mKJ 

V / 

LC 


(0.454m) 


/ 


9.476  — 
K 


= 0.0011 


C 

— osjl  = 0.441 


NTU=0.368 


(4.3.3. 1.3) 

(4.3.3.1.4) 

(4.3.3. 1.5) 


These  parameters  lie  below  the  range  of  data  given  by  Ranganayakulu.  Using 
extrapolation,  it  was  determined  that  the  conduction  effect  factor  would  be  approximately 
0.0005.  This  means  that  this  particular  tube  in  the  heat  exchanger  would  see  a loss  in 
capacity  of  one  twentieth  of  one  percent  due  to  axial  heat  conduction.  When  this  capacity 
degradation  is  summed  over  all  of  the  tubes  in  the  entire  heat  exchanger  where  this  effect 
occurs,  the  capacity  reduction  totals  0.13  W (0.45  Btu/h),  which  is  insignificant  when 
compared  to  the  predicted  performance  of  the  evaporator  being  49  000  W (16  800  Btu/h). 

It  should  be  noted  that  the  effectiveness  of  tube  15  is  0.29.  Hence,  our  result  agrees  with 
the  general  statements  by  Kays  and  London  and  that  of  Ranganayakulu  et.  al.  that  the 
longitudinal  heat  transfer  has  an  insignificant  effect  for  heat  exchangers  with  an 
effectiveness  less  than  0.8.  The  negligible  impact  of  the  longitudinal  tube  conduction  on 
evaporator  performance  permits  neglecting  this  heat  transfer  in  modeling  of  a finned-tube 
heat  exchanger.  It  may  be  further  inferred  that  the  same  conclusion  can  be  made  for  the 
R407C  zeotropic  mixture.  Although  a 7 °C  (12  °F)  glide  associated  with  R407C  phase 
change  produces  a temperature  difference  promoting  longitudinal  heat  transfer,  in  the 


70 


analyzed  evaporator  this  glide  would  be  distributed  over  10  m (33  ft)  of  tube  passes  and 
would  result  in  a small  longitudinal  temperature  gradient. 

4.3.S.2  Tube-to-Tube  Heat  Transfer  via  Fins 

If  we  recognize  that  longitudinal  tube  heat  conduction  has  a negligible  impact,  then  the 
difference  in  capacity  degradation  between  COIL-E  and  COIL-EC  at  16.7  °C  (30  °F) 
superheat  must  be  due  to  longitudinal  fin  conduction.  In  COIL-EC,  the  continuous  cuts  in 
the  fins  separating  different  depth  rows  prevent  heat  transfer  between  different  depth 
rows.  However,  the  fins  join  the  adjacent  tubes  in  the  same  depth  row,  and  some  heat 
transfer  between  them  occurs.  This  is  why  COIL-EC  still  experiences  a decline  in 
capacity,  but  not  as  much  as  COIL-E. 

Our  literature  review  located  two  publications  that  shed  some  light  on  the  longitudinal  fin 
conduction  phenomenon.  Heun  and  Crawford  (1994)  performed  analytical  study  of  the 
effects  of  longitudinal  fin  conduction  on  multipass  cross-counterflow  single-depth-row 
heat  exchanger.  They  considered  the  fins  to  have  one-dimensional  temperature 
distributions  and  solve  them  using  a system  of  non-dimensional  differentials  equations. 
Their  results  showed  that  longitudinal  fin  conduction  always  degrades  heat  exchanger 
performance  and  this  effect  is  stronger  for  a low  normalized  fin  resistance  and  large 
values  of  the  ratio  of  air-side  conductance  to  air  heat  capacity  rate. 

Romero-Mender  et  al.  (1997)  also  studied  tube-to-tube  heat  transfer  in  a single-row 
finned-tube  heat  exchanger.  They  assumed  the  fins  to  be  continuously  and  uniformly 
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distributed  along  the  length  of  each  tube.  With  his  continuum  assumption,  they  solved  a 
system  of  ordinary  differential  equations  for  steady-state  refrigerant  and  tube-wall 
temperatures.  They  identified  four  non-dimensional  groups  that  effected  the  degradation 
of  evaporator  capacity.  These  group  are:  (1)  the  ratio  of  the  thermal  conductance  for 
convective  heat  transfer  between  the  refrigerant  and  the  wall  to  the  product  of 
refrigerants  heat  capacity  and  mass  flow  rate,  (2)  the  ratio  of  the  thermal  conductance  for 
external  heat  transfer  from  the  unfinned  portion  of  the  tube  to  the  internal  thermal 
conductance,  (3)  the  ratio  of  the  thermal  conductance  for  convection  from  the  fin  to  the 
thermal  conductance  for  conduction  along  it,  and  (4)  the  ratio  of  the  thermal  conductance 
for  heat  conduction  along  the  insulated  fin  to  the  thermal  conductance  between  the 
refrigerant  and  the  wall.  Their  study  also  indicated  the  number  of  tubes  to  be  an 
influencing  factor.  The  study  by  Romero-Mender  et  al.  indicates  that  tube-to-tube  heat 
transfer  always  degrades  capacity  and  that  the  influencing  parameters  they  identified 
have  a non-linear  impact  on  capacity  degradation  over  the  wide  range  of  values  studied. 
For  some  parametric  values  they  found  the  degradation  in  a single-row  heat  exchanger  to 
be  as  high  as  20  %. 

Our  literature  review  have  not  located  a publication  that  would  discuss  longitudinal  fin 
conduction  in  a multi-depth-row  evaporator,  but  it  may  be  safely  expected  that  capacity 
degradation  would  be  higher  than  that  for  a single-depth-row  coil.  The  literature  review 
identified  a very  recent  paper  which  presents  three  simulation  models  for  finned-tube 
single-phase  dehumidifying  heat  exchangers,  the  most  advanced  of  which  accounts  for 
tube-to-tube  heat  transfer  (Oliet  et  al.  (2002)).  This  model,  referred  to  as 
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“advancedCHESS”,  is  based  on  a finite  volume  approach  with  discretization  of  the  heat 
exchanger  domain  into  a set  of  control  volumes  as  fin-and-tube  elements  where  both 
local  thermophysical  properties  and  empirical  coefficients  are  computed.  While 
“advancedCHESS”  appears  to  be  a research  model,  two  other  models  are  more  practical 
for  production  simulations.  The  two  less  advanced  models,  called  “basicCHESS”  and 
“quickCHESS”,  do  not  consider  tube-to-tube  heat  transfer. 

While  the  above  publications  are  very  interesting  and  valuable,  they  do  not  offer  a 
methodology  for  accounting  for  tube-to-tube  heat  transfer  in  a tube-by-tube  simulation 
model.  The  number  of  influencing  parameters  identified  for  a dry  fin  by  Romero- 
Mender  et  al.  (1997)  suggests  that  a fully  fundamental  approach  will  be  difficult  to 
implement  into  a tube-by-tube  evaporator  model,  which  uses  the  adiabatic  fin  tip 
assumption  and  considers  an  individual  tube  as  a separate  entity  for  heat  transfer 
calculations.  Our  attempts  to  apply  a few  algorithms  derived  from  their  paper,  however, 
were  unsuccesful  because  we  were  unable  to  resolve  the  “clashes”  between  the 
fundamental  algorithms  and  the  current  simulation  scheme  used  in  EVAP5.  It  appears 
that  merging  a fundamental  scheme  into  EVAP5  amounts  to  a separate  project  that 
should  be  dedicated  to  this  task. 

To  reach  the  objectives  of  the  project  within  a stipulated  effort  and  time,  a practical 
scheme  was  developed,  which  uses  the  temperature  difference  between  neighboring  tubes 
as  the  driving  force  for  heat  transfer.  This  scheme  approaches  the  tube-to-tube  heat 
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transfer  problem  in  a similar  way  Sheffield  (1988)  studied  fin  collar-tube  heat  transfer 
resistance  as  shown  in  Figure  4.3.3.2.I. 


Figure  4. 3. 3.2.1 : Schematic  graph  for  longitudinal  fin  conduction  between  two  adjacent 
tubes 


To  determine  the  heat  transferred,  the  Fourier  Law  of  Conduction  is  applied.  The  effects 
of  the  available  width  and  configuration  of  the  conducting  material  (fin)  are  represented 
by  a “shape  factor”  S used  in  equation  4.9: 


( Qfin)i,j  ~ 


W-tf 


Kf 


(Twj  -T*.i)=Stj-  K, (TWJ  - TWJ ) 


(4.3. 3.2.1) 


To  show  the  impact  of  fin  conduction,  Lee  and  Domanski  used  this  scheme  considering 
up  to  six  immediate  neighboring  tubes.  For  example,  for  the  circuitry  used  in  the 
evaporators  tested  under  this  project  and  shown  in  Figure  4. 3. 3. 2. 2 the  intermediate 


neighbors  for  tube  25  are  tubes  7,  8,  26,  44,  43,  and  24. 
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Figure  4. 3. 3. 2. 2:  Schematic  of  refrigerant  circuitry  for  COIL-W,  COIL-E  and  COIL-EC 
in  cross-counter  flow  configuration  (inlet  tubes:  42,  48,  54;  outlet 
tubes:  1,7,  13) 


Extensive  experimenting  with  this  scheme  for  the  coils  tested  under  this  project  indicated 
that  it  was  important  to  add  additional  neighbors  to  the  group  of  immediate  neighbors 
considered  so  far.  This  need  demonstrated  itself  not  only  in  predicted  capacity  values 
but  also  in  simulation  runs,  which  did  not  yield  gradually  changing  predictions  at  small 
changes  in  imposed  refrigerant  superheat  at  the  outlet  tubes.  Based  on  these 
observations,  depending  on  location  up  to  six  second-order  neighbors  were  added.  These 
are  other  tubes  in  the  coil  assembly  that  a given  tube  “can  see”.  Tube  25  has  four 
second-order  neighbors;  they  are  tubes  6,  9,  45,  and  42.  For  tube  9,  the  immediate 
neighbors  are  tube  8,  10,  27,  and  26,  and  the  second-order  neighbors  are  tubes  25,  45, 
and  28.  In  a three-depth-row  coil,  the  maximum  number  of  second-order  neighbors  is 
four.  A five-depth-row  coil  is  need  for  a tube  located  in  the  middle  depth  row  to  have  all 
six  second-order  neighbors. 


The  value  of  the  shape  factor  depends  on  a fin  design.  For  flat  and  wavy  fins  the  fin 
material  is  continuous.  Lanced  fins,  however,  have  numerous  cuts,  which  reduce  the  fin 
cross-section  area  that  is  available  for  heat  transfer.  Hence,  the  shape  factor  for  flat  and 
wavy  fins  should  be  expected  to  have  higher  values  than  for  lanced  or  louver  fins.  Since 
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we  are  not  aware  of  any  publication  that  quantifies  the  fin  shape  factor,  the  values  for 
COIL-W  and  COIL-E  were  assigned  based  on  their  respective  results  for  test  10.  Once 
these  values  for  shape  factors  were  assigned,  they  were  left  unchanged  for  the  remaining 
simulations.  COIL-EC  used  the  same  shape  factor  as  COIL-E,  however,  any  tube  in 
COIL-EC  could  have  only  two  neighbors,  the  closest  two  tubes  located  in  the  same  depth 
row. 


Figure  4. 3. 3. 2. 3 shows  tested  and  simulated  capacities  for  COIL-E  at  conditions  of  test  1, 
2,  9,  10,  13,  and  14.  For  the  tests  at  5.6  °C  (10  °F)  superheat  (tests  1,  9 and  13),  the 
model  predicted  measured  capacities  within  5.1  %.  For  the  tests  with  16.7  °C  (30  °F) 
superheat,  the  differences  in  between  tested  and  predicted  and  tested  capacities  were 
- 6.7  %,  3.1  %,  and  - 2.9  %.  Without  accounting  for  tube-to-tube  heat  transfer,  EVAP5 
would  overpredict  the  capacities  at  16.7  °C  (30  °F)  superheat  by  approximately  20  %. 
Section  4.4  presents  validation  results  for  all  three  evaporators. 


Figure  4. 3. 3. 2. 3:  COIL-E  measured  and  simulated  capacities  for  tests  with  5.6  °C  (10  °F) 
and  16.7  °C  (30  °F)  superheats  (tests  1,  2,  9,  10,  13,  14) 
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4.4  Validation  and  Simulations  with  EVAP5 


4.4. 1 Validation  of  EVAP5 

The  majority  of  laboratory  test  data  measured  in  this  study  are  for  a cross-counter  flow 
configuration  with  non-disturbed  (uniform)  air  velocity  profile  at  the  coil  inlet.  These 
measurements  for  COIL-W,  COIL-E,  and  COIL-EC  were  used  to  validate  EVAP5  and 
explore  the  impact  of  tube-to-tube  heat  transfer.  The  imposed  variety  of  refrigerant 
superheats  at  individual  exit  tubes  combined  with  a wide  range  of  air  flow  provided  a 
unique  set  of  challenging  test  data  for  validating  any  evaporator  model. 

The  refrigerant  circuitry  for  the  tested  evaporators  has  already  been  presented  in  Figure 
4. 3. 3. 2. 2 as  it  is  displayed  by  the  EVAP-COND  interface.  Also  copying  from  the 
respective  window  of  EVAP-COND,  Figure  4.4. 1.1  shows  the  key  design  parameters  of 
COIL-W.  Except  for  a different  fin  design,  these  parameters  were  the  same  for  COIL-E 
and  COIL-EC. 

At  the  outset  of  simulations  for  each  coil,  EVAP5  was  “tuned”  to  predict  the  performance 
of  a given  evaporator  at  the  conditions  of  test  9.  This  was  accomplished  by  inputting 
“corrections  parameters”  for  the  refrigerant  heat  transfer  coefficient,  refrigerant  pressure 
drop,  and  air-side  heat  transfer  coefficient.  (Section  4.3.2  discusses  the  reasons  for  using 
these  parameters  in  the  context  of  prediction  discrepancies  between  different  air-side 
correlations).  Figure  4.4. 1.2  presents  the  correction  parameters  for  COIL-W  and  COIL-E 
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as  they  were  input  into  the  EVAP-COND  window.  The  input  for  COIL-EC  was  different 
by  the  value  for  the  air-side  heat  transfer  coefficient,  which  was  0.62  instead  of  0.65.  The 
1 .6  value  for  the  refrigerant  pressure  drop  parameter  accounts  for  the  impact  of  lubricant, 
which  can  be  responsible  for  35  % pressure  drop  underprediction.  Since  the  parameter 
for  refrigerant  heat  transfer  coefficient  was  set  to  1.0,  the  0.65  or  0.62  value  for  the  air- 
side  heat  transfer  coefficient  accommodates  the  heat  transfer  adjustment  on  the  air  and 
refrigerant  sides.  These  correction  parameters  were  used  in  all  simulations  for  the 
respective  coils. 


Figure  4.4. 1.1:  Design  parameters  for  COIL-W 
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Figure  4. 4. 1.2:  Correction  parameters  for  COIL-W  and  COIL-E 

Tables  4. 4. 1.1,  4. 4. 1.2,  and  4. 4. 1.3  show  tested  and  simulated  total  capacities  for  COIL- 
W,  COIL-E,  and  COIL-EC,  respectively,  in  the  test  matrix  format.  Tables  4.4.1.4a, 
4.4.1.4b,  4.4.1.5a,  4.4.1.5b,  4.4.1.6a,  and  4.4.1.6b  present  total  and  sensible  capacities, 
sensible  heat  ratios,  and  differences  between  simulated  and  measured  results.  It  is 
convenient  to  screen  the  accuracy  of  capacity  predictions  by  reviewing  Figures  4. 4. 1.3 
and  4. 4. 1.4.  Figure  4. 4. 1.3  shows  that  the  maximum  error  in  capacity  prediction  for  wet 
coil  tests  was  5.5  % for  all  air  velocities  and  refrigerant  superheat  scenarios.  For  dry  coil 
tests  shown  in  Figure  4. 4. 1.4,  EVAP5  predicted  the  three  capacities  at  uniform  5.6  °C 
(10  °F)  superheat  within  5.0  %.  (They  are  represented  by  the  first  bar  for  each  coil).  For 
the  remaining  six  tests  at  different  superheat  scenarios,  the  only  capacity  predicted  within 
5.0  % was  for  COIL-EC  (represented  by  the  last  bar  on  the  right  hand  side).  Capacities 
for  COIL-W  and  COIL-E  capacities  were  unpredicted  by  as  much  as  20.0  %.  The 
inability  of  EVAP5  to  account  accurately  for  longitudinal  fin  conduction  for  dry  coil  tests 
can  be  explained  by  the  fact  that  the  used  algorithm  considers  only  temperature 
differences  between  neighboring  tubes  (the  first  order  effect)  and  neglects  variations  in 
air-side  heat  transfer.  This  conclusion  agrees  with  the  observation  by  Romero-Mendez  et 
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al.  (1997)  who  indicated  thermal  conductance  for  convection  from  the  fin  as  one  of  the 
parameters  affecting  tube-to-tube  heat  transfer.  This  and  other  effects  identified  by 
Romero-Mendez  et  al.  (1997)  should  receive  detailed  attention  in  a future  study  dedicated 
to  this  challenging  modeling  issue. 
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Table  4.4. 1.1:  Measured  and  Simulated  Capacities  for  COIL-W  in  Cross-Counter  Flow  Configuration 


oo 


* Superheat  to  be  controlled  such  that  the  desired  overall  level  of  superheat  is  obtained 

1)  SI  units  of  mVkWh  multiplied  by  capacity  (Q)  in  kW  to  determine  airflow  (IP  units  of  cfm/ton  multiplied  by  capacity  (Q)  in  tons). 


Table  4.4.1 .2:  Measured  and  Simulated  Capacities  for  COIL-E  in  Cross-Counter  Flow  Configuration 


Cn| 
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Superheat  to  he  controlled  such  that  the  desired  overall  level  of  superheat  is  obtained 

SI  units  of  m3/kWh  multiplied  by  capacity  (Q)  in  kW  to  determine  airflow  (IP  units  of  cfm/ton  multiplied  by  capacity  (Q)  in  tons) 


Table  4. 4. 1.3:  Measured  and  Simulated  Capacities  for  COIL-EC  in  Cross-Counter  Flow  Configuration 


m 

oo 


1)  SI  units  of  m3/kWh  multiplied  by  capacity  (Q)  in  kW  to  determine  airflow  (IP  units  of  cfm/ton  multiplied  by  capacity  (Q)  in  tons). 


Table  4,4. 1 .4a:  EVAP5  Validations  for  COIL-W  (evaporator  with  wavy  fins),  SI  Units 
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Table  4.4.1.4b:  EVAP5  Validations  for  COIL-W  (evaporator  with  wavy  fins),  IP  Units 
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Table  4.4.1.5a:  EVAP5  Validations  for  COIL-E  (evaporator  with  lanced  fins),  SI  Units 
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Table  4.4,1.5b:  EVAP5  Validations  for  COIL-E  (evaporator  with  lanced  fins),  IP  Units 
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Table  4.4. 1 ,6a:  EVAP5  Validation  for  COIL-EC  (evaporator  with  lanced  fins,  cut  between  tube  depth  rows),  SI  Units 
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Figure  4.4.1 .3:  Difference  between  simulated  and  measured  capacities  for  all  wet  coil  tests  for 
COIL-W,  COIL-E,  and  COIL-EC 


Figure  4.4. 1.4  Difference  between  simulated  and  measured  capacities  for  all  dry  coil  tests  for 
COIL-W,  COIL-E,  and  COIL-EC 
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Tube-to-tube  heat  transfer  demonstrates  itself  in  temperatures  that  can  be  measured  on  coils 
return  bends,  as  it  was  shown  for  COIL-W  in  Figures  3.4.4. 1 and  3 .4.4.2.  Figure  4.4. 1.5  shows 
similar  information  (refrigerant  temperature  at  tube  exits)  as  it  is  displayed  by  EVAP-COND  for 
the  same  tests.  For  test  9 with  even  refrigerant  superheat,  refrigerant  temperatures  are  similar  for 
each  circuit;  refrigerant  temperatures  reflect  drop  in  refrigerant  pressure  until  the  last  two  tubes 
in  each  circuit  (2  and  1,  8 and  7,  and  14  and  13)  in  which  the  refrigerant  is  superheated.  For  test 
12,  the  first  circuit  is  in  two-phase  flow  until  the  exit  tube  1,  while  the  refrigerant  leaving  two 
other  exit  tubes  (7  and  13)  is  highly  superheated.  Tubes  7 and  8 experience  a drop  in 
temperature  compared  to  tube  9 because  of  their  vicinity  to  the  left-hand  side  circuit  with  two- 
phase,  low-temperature  refrigerant.  Tubes  13  and  14  also  experience  temperature  drop, 
however,  it  is  small  because  the  adjacent  tubes  are  also  superheated.  This  simulation  results 
agree  in  principle  with  the  measured  return  bend  temperature  of  Figures  3.4.4. 1 and  3.4.4.2. 
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Tube  numbering 


A 


air  flow 


Refrigerant  exit  temperatures  (°F)  for  COIL-W,  test  9 


Refrigerant  exit  temperatures  (°F)  for  COIL-W,  test  12 


Figure  4.4. 1 .5:  Tube  numbering  and  refrigerant  exit  temperatures  for  individual  tubes  for  test  9 
and  test  12  for  COIL-W,  (inlet  tubes:  42,  48,  54;  outlet  tubes:  1,7,  13) 


Figure  4.4. 1.6  presents  refrigerant  exit  qualities  for  individual  tubes  for  COIL-E  and  COIL-EC 
for  test  10  (16.7  °C  (30  °F)  even  superheat).  The  figure  shows  that  refrigerant  reaches  superheat 
(quality  1)  two  tubes  earlier  in  COIL-E  than  in  COIL-EC  in  each  refrigerant  circuit.  This  is  a 
result  of  heat  transfer  between  different  tube  depth  rows  that  was  allowed  in  COIL-E  and  was 
inhibited  in  COIL-EC.  Corresponding  tube  temperatures  predicted  by  EVAP5  agree  with  those 
measured  during  laboratory  tests. 
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Refrigerant  exit  qualities  (fraction)  for  COIL-E,  test  10 


Refrigerant  exit  qualities  (fraction)  for  COIL-EC,  test  10 


Figure  4.4. 1.6:  Refrigerant  exit  qualities  for  individual  tubes  for  COIL-E  and  COIL-EC  test  10 
(16.7  °C  (30  °F)  even  superheat) 


The  last  two  columns  in  Tables  4.4. 1.4,  4.4. 1.5,  and  4.4. 1.6  compare  simulated  capacities  that 
were  obtained  with  and  without  accounting  for  tube-to-tube  heat  transfer.  For  the  tests  with  a 
uniform  superheat  of  5.6  °C  (10  °F),  the  difference  in  capacities  is  not  greater  than  2.7  % for  any 
of  the  three  coils.  For  the  tests  with  uniform  superheat  of  16.7  °C  (30  °F),  the  capacities  differ 
by  63.1  % and  56.3  % for  COIL-W,  33.6  %,  18.8  %,  19.2  % for  COIL-E,  and  1.5  %,  2.5  %, 
2.5  %,  and  4.9  % for  COIL-EC.  These  results  demonstrate  the  impact  of  fin  design  on  tube-to- 
tube  heat  transfer. 


Thevalidation  of  EVAP5  used  the  full  set  of  COIL-W,  COIL-E,  and  COIL-EC  measurements  in 
cross-counter  flow  configuration,  which  constituted  the  majority  of  the  measurements  taken  in 
this  study.  The  validation  effort  was  not  extended  to  the  eight  data  points  taken  in  cross-parallel 
flow  configuration  because  six  of  these  tests  that  involved  16.7  °C  (30  °F)  superheat  resulted  in 
severe  pinching  within  less  than  1 °C  (1.8  °F)  in  at  least  one  of  the  circuits.  Such  a close 
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approach  of  refrigerant  and  air  causes  a profound  convergence  problem  for  a tube-by-tube  model 
in  which  air  temperatures  upstream  of  each  tube  have  to  be  iterated  around  a target  value. 
Furthermore,  for  evaluating  the  potential  reduction  in  heat  exchanger  volume,  only  capacity 
predictions  at  test  9 with  5.6  °C  (10  °F)  superheat  were  needed,  and  these  were  attainable  with 
EVAP5. 


Test  9 measured  capacities  for  COIL-W  and  COIL-E  were  4729  W (16146  Btu/h)  and  4549  W 
(15522  Btu/h),  while  EVAP5  predictions  were  4796  W (16366  Btu/h)  and  5482  W (18705 
Btu/h),  respectively.  This  is  a very  good  prediction  for  COTL-W,  within  1.3  %,  while  the 
discrepancy  for  COIL-E  is  20.5  %.  It  should  be  noted  that  a higher  capacity  should  be  expected 
for  COIL-E  than  for  COIL-W,  as  it  was  predicted  by  EVAP5  and  always  was  obtained  from 
laboratory  measurements  except  this  time.  It  is  possible  that  some  condensate  holdup  might 
have  influenced  the  measured  capacity  for  COIL-E.  With  this,  it  was  concluded  that  EVAP5 
properly  simulated  coils  in  a cross-parallel  flow  set  up.  Consequently,  COIL-W  was  applied  in  a 
later  section  to  examine  potential  savings  in  evaporator  core  volume  for  the  cross-parallel 
configuration. 

4.4.2  Possible  Savings  in  Heat  Transfer  Area  Due  to  Optimized  Superheat 

4.4.2. 1 Cross-Counter  Flow  Configuration  with  Uniform  Air  Flow  Distribution 
Considering  similar  performance  degradations  for  different  refrigerant  superheat  scenarios, 
possible  savings  in  heat  exchanger  material  are  demonstrated  using  test  12  of  COIL-W  as  an 
example.  In  our  simulations,  it  was  assumed  that  smart  refrigerant  distributors  would  optimize 
refrigerant  distribution  so  the  evaporator  obtains  maximum  capacity.  In  these  tests  with  a 
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uniform  air  velocity  profile,  the  evaporators  reached  maximum  capacity  when  the  refrigerant 
split  between  the  three  circuits  resulted  in  uniform  superheat  at  the  individual  outlet  tubes. 

For  these  simulations,  five  alternative  coils  were  coded  with  a smaller  number  of  tubes  than 
COIL-W.  All  simulation  runs  had  the  same  inlet  air  condition,  refrigerant  inlet  quality,  and 
refrigerant  outlet  pressure  and  superheat  at  the  evaporator  exit.  Also,  inlet  air  velocity  was  the 
same  for  each  coil  as  for  COIL-W.  A coil  with  a smaller  face  area  had  a lower  volumetric  flow 
rate  than  COIL-W,  proportional  to  the  percentage  that  its  face  area  was  reduced. 

Figure  4.4. 2. 1.1  shows  coil  designs  and  simulation  results.  Four  out  of  five  alternative  coil 
designs  offered  both  savings  in  the  heat  exchanger  core  volume  and  an  increase  in  coil  capacity. 
The  coils  with  a lower  volumetric  flow  of  air  would  also  provide  savings  in  fan  power. 
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Figure  4.4. 2. 1.1:  Simulation  results  for  alternative  coil  designs  to  COIL-W  in  cross-counter  flow 
configuration.  Performance  is  compared  to  COIL-W  test  12  with  simulated 
capacity  of  13225  Btu/h. 


4. 4. 2. 2 Cross-Parallel  Flow  Configuration  with  Uniform  Air  Flow  Distribution 
Simulations  were  also  performed  to  demonstrate  possible  savings  in  coil  material  (core  volume 
of  the  heat  exchanger)  for  COIL-W  in  cross-parallel  configuration.  Test  9 (W020304a)  was  used 
as  a reference.  Alternative  coil  designs  with  two  depth  rows  were  only  examined  because,  in  the 
cross-parallel  configuration,  more  depth  rows  are  not  beneficial  due  to  pinching. 
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All  simulations  were  run  using  test  9 operating  conditions,  including  6 °C  (10.8  °F)  superheat, 
with  the  difference  that  the  volumetric  flow  of  air  was  adjusted  so  that  each  coil  had  the  same 
inlet  air  velocity  as  COIL-W.  This  means  that  a coil  with  a smaller  face  area  had  a lower 
volumetric  flow  rate  than  COIL-W  by  the  same  percentage  its  face  area  was  reduced. 

Figure  4.4. 2. 2.1  shows  coil  designs  and  simulation  results.  Each  of  the  four  presented  two-depth 
row  designs  offered  improved  capacity  and  savings  in  coil  core  volume.  The  coils  with  a lower 
volumetric  flow  of  air  would  also  provide  savings  in  fan  power.  The  smallest  coil  with  12x2 
tube  arrangement  matched  the  capacity  of  test  12  with  a 33.3  % savings  in  coil  material. 
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Figure  4.4.2.2.1:  Simulation  results  for  alternative  coil  designs  to  COIL-W  in  cross-parallel 
configuration.  Comparisons  are  to  COIL-W,  test  12  with  tested  capacity  of 


7311  Btu/h. 


4. 4. 2. 3 Cross-Counter  Flow  Arrangement  with  Non-Uniform  Air  Flow  Distribution 
The  tests  performed  in  the  lab  with  non-uniform  air  distributions  resulted  in  complicated  velocity 
profiles  that  currently  cannot  be  reproduced  in  EVAP5  simulations.  For  this  reason,  the 
simulations  were  performed  with  one-dimensional,  non-uniform  velocity  profiles  that  were 
independent  of  the  tests  performed  in  the  lab  and  represented  a different  application  case 
scenario.  For  these  simulations,  COIL-W  test  9 with  a uniform  air  distribution  was  selected  as  a 
reference  test,  and  additional  simulations  were  performed  for  two-step  velocity  profiles  where 
the  top  (left)  to  bottom  (right)  velocity  ratios  were  1:1.5,  1 :2.0,  1 :2.5,  1 :3,  1:3.5,  and  1 :4. 
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Two  runs  were  performed  for  each  velocity  ratio:  the  first  - with  a uniform  refrigerant 
distribution,  and  the  second  - with  a refrigerant  distribution  optimized  to  obtain  maximum 
capacity.  During  all  these  tests,  the  external  run  parameters  (refrigerant  inlet  state,  exit  pressure 
and  superheat,  air  flow  rate,  etc.)  were  the  same.  Figure  4.4. 2. 3.1  presents  a velocity  profile 
representation  for  the  1 :3  ratio  as  it  was  input  into  EVAP-COND.  Since  our  velocity  profiles 
have  a near  step  change,  they  represent  a more  radical  departure  from  uniformity  than  the 
profiles  obtained  in  the  laboratory. 

Table  4. 4. 2. 3.1  summarizes  simulation  results,  and  Figure  4.4. 2. 3. 2 presents  simulated  capacities 
as  referenced  to  the  capacity  of  test  9.  The  table  and  figure  show  that  the  capacity  degrades 
linearly  with  degradation  of  the  air  velocity.  For  the  1:4  air  velocity  ratio  and  uniform 
refrigerant  distribution,  the  obtained  capacity  was  only  63  % of  the  reference  test  9 value. 
However,  with  optimized  refrigerant  distribution,  as  is  the  purpose  of  smart  distributors,  the 
obtained  capacity  was  within  7 % of  the  reference  capacity. 
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Figure  4.4. 2. 3. 1 : Air  velocity  profile  representation  for  the  1 :3  top-to-bottom  velocity  ratio 
(during  tests  the  coil  was  positioned  vertically,  turned  clockwise  by  90  °) 


Table  4.4. 2. 3.1 : Simulated  Capacities  and  Refrigerant  Distributions  for  Non-Uniform  Inlet  Air 


Velocity  Profile 


Air  velocity  ratio  (top  to  bottom) 

1 : 1 

1 : 1.5 

1 : 2 

1 : 2.5 

1 : 3 

1 : 3.5 

1 : 4 

Capacity 

(watt) 

Uniform  ref  distribution 

7044 

6748 

6194 

5710 

5199 

4886 

4465 

Optimized  ref  distribution 

7044 

7001 

6914 

6849 

6762 

6662 

6582 

Capacity 

(Btu/h) 

Uniform  ref.  distribution 

22127 

21197 

19456 

17935 

16331 

15347 

14024 

Optimized  ref.  distribution 

22127 

21997 

21718 

21515 

21240 

20928 

20675 

Optimized  refrig, 
distribution 
(fraction) 

top  (left)  circuit 

0.33 

0.295 

0.265 

0.240 

0.222 

0.210 

0.195 

middle  circuit 

0.33 

0.333 

0.340 

0.340 

0.342 

0.345 

0.350 

bottom  (right)  circuit 

0.33 

0.372 

0.395 

0.420 

0.436 

0.445 

0.455 
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1.100 


Air  velocity  ratio  (bottom/top) 


Figure  4.4. 2. 3. 2:  COIL-W  capacities  at  different  air  velocity  ratios  referenced  to  capacity  at  test 
9 in  cross-counter  flow  configuration 


To  assess  the  savings  in  the  heat  exchanger  material  due  to  optimized  control  of  refrigerant 
superheat,  two  evaporators  with  reduced  number  of  tubes  were  coded,  shown  in  Figure  4.4.2. 3. 3, 
and  simulations  were  performed  for  1:2.5  and  1:4  air  velocity  profiles.  The  two-depth  row 
evaporator  had  the  same  face  area  as  COIL-W  and  was  simulated  with  the  same  volumetric  flow 
rate  of  air.  For  the  three-depth-row,  the  volumetric  flow  rate  was  reduced  by  16.7  %,  which 
corresponds  to  the  reduction  of  the  coil  face  area  in  relation  to  that  of  COIL-W.  The  results 
presented  in  Table  4. 3. 3. 1.1  show  that  the  benefit  of  optimizing  refrigerant  distribution  increases 
with  the  level  of  non-uniformity  in  the  air  velocity  profile.  For  the  1:4  air  velocity  ratio, 
optimizing  refrigerant  distribution  allows  a reduction  in  coil  volume  of  33.3  %.  For  the  1 :2.5  air 
velocity  ratio,  the  use  of  15x3  coil  with  a slightly  increased  volumetric  flow  rate  could  produce  a 
16.7  % reduction  in  the  coil  volume. 


100 


Coil  15x3 


Coil  18x2 


Figure  4.4. 2. 3. 3:  Two  evaporators  with  a reduced  number  of  tubes 


Table  4.4. 2. 3. 2:  Savings  in  Coil  Volume  Relative  to  COIL-W  due  to  Optimized  Refrigerant  for 
1 :2.5  and  1 :4  Air  Velocity  Ratios 


Coil 

Refrigerant 

distribution 

Air  velocity 
ratio 

(top  to  bottom) 

Air 

volumetric 
flow  rate 
mVmin 

Capacity 

watt 

Air 

volumetric 
flow  rate 
ft3/min 

Capacity 

Btu/h 

Savings 
in  coil 
volume 

% 

COIL-W 

uniform 

1 : 1 

20.7 

6485 

733 

22127 

0 

COIL-W 

optimized 

1 : 2.5 

20.7 

6305 

733 

21515 

0 

COIL-W 

uniform 

1 : 2.5 

20.7 

5256 

733 

17935 

0 

18x2 

optimized 

1 : 2.5 

20.7 

4558 

733 

15553 

33.3 

15x3 

optimized 

1 : 2.5 

17.3 

4890 

611 

16685 

16.7 

COIL-W 

uniform 

1 : 1 

20.7 

6485 

733 

22127 

0 

COIL-W 

optimized 

1 : 4 

20.7 

6059 

733 

20675 

0 

COIL-W 

uniform 

1 : 4 

20.7 

4110 

733 

14024 

0 

18x2 

optimized 

1 : 4 

20.7 

4335 

733 

14790 

33.3 

15x3 

optimized 

1 : 4 

17.3 

4474 

611 

15266 

16.7 
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5 CONCLUSIONS 


This  collection  of  experimental  data  for  the  three  evaporators  has  revealed  interesting  results 
related  to  non-uniform  refrigerant  distribution  and  conduction  between  tubes  through  the  fins. 
With  cross-counter  refrigerant  flow,  uniform  airflow,  and  exit  manifold  superheat  fixed  at  5.6  °C 

(10.0  °F),  the  wavy  fin  and  wavy-lanced  fin  evaporator’s  capacity  dropped  by  as  much  as  41  % 
and  32  %,  respectively,  as  the  superheat  was  allowed  to  vary  between  the  circuits.  Control  of 
superheat  was  shown  to  be  even  more  important  during  cross-parallel  refrigerant  flow  due  to  the 
rapid  pinching  of  the  refrigerant  and  air  temperatures.  For  the  wavy  and  lanced  finned 
evaporators  in  cross-parallel  flow,  capacity  dropped  by  85  % and  78  % as  superheat  changed 
from  5.6  °C  (10.0  °F)  to  16.7  °C  (30.0  °F). 

As  the  coil’s  faces  were  blocked  to  produce  a non-uniform  airflow,  control  of  superheat  was 
shown  to  restore  capacity  if  the  volumetric  flow  of  air  was  unchanged.  The  tests  showed  that 
when  airflow  rate  was  held  constant,  the  losses  in  capacity  due  to  non-uniform  airflow  could  be 
recovered  to  within  2 % of  the  original  uniform  airflow  capacity  by  controlling  superheat.  The 
more  non-uniform  the  airflow  over  the  coil,  the  greater  was  the  benefit  of  controlling  superheat. 
For  the  lanced  fin  coil,  as  the  airflow  ratio  between  the  top  half  and  lower  half  of  the  coil  varied 
from  1:1.26  to  1:2.59,  superheat  control  improved  capacity  by  1.4  % and  4.6  %,  respectively. 

In  parallel  with  the  experimental  effort,  the  NIST  evaporator  model  EVAP5  was  upgraded  to 
control  refrigerant  distribution  and  account  for  tube-to-tube  heat  transfer.  The  model  was 
validated  with  the  experimental  results  and  then  used  to  determine  the  possible  savings  in 
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evaporator  core  volume  if  refrigerant  distribution  was  controlled  by  a smart  distributor.  In 
extreme  cases,  the  savings  in  core  volume  could  be  as  much  as  40  %. 

A combination  of  results  obtained  from  laboratory  testing  and  simulations  indicated  the  influence 
of  tube-to-tube  heat  transfer  on  capacity  degradation.  The  impact  of  tube-to-tube  heat  transfer 
was  negligible  in  tests  with  a uniform  5.6  °C  (10  °F)  superheat,  but  it  was  significant  in  tests 
involving  16.7  °C  (30  °F)  superheat.  Between  two  possible  conduction  mechanisms  by  which 
such  heat  transfer  may  occur,  longitudinal  fin  conduction  was  chiefly  responsible  for  degraded 
performance  while  longitudinal  tube  conduction  had  insignificant  effect.  The  upgraded  version 
of  the  EVAP5  evaporator  model,  which  accounts  for  tube-to-tube  heat  transfer  based  on  tube 
temperatures,  was  able  to  predict  key  return  bend  temperatures  that  indicated  the  occurrence  of 
tube-to-tube  heat  transfer.  However,  the  study  also  confirmed  that  longitudinal  heat  conduction 
is  affected  by  the  fin  design,  air-side  heat  transfer  coefficient,  and  moisture  removal  process. 
Consequently,  a more  detailed  modeling  scheme  needs  to  be  developed  to  capture  other  effects 
influencing  tube-to-tube  heat  transfer.  Such  a study  would  not  only  improve  the  modeling  of 
evaporators  but  also  of  condensers  and  of  gas  coolers,  where  internal  heat  transfer  may  be  even 
more  pronounced. 
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APPENDIX  A.  SUMMARY  OF  TEST  RESULTS 


A.  1 Wavy  fin  evaporator  in  cross-counter  flow 


Table  A.  1 . 1 : Wavy  Fin  Evaporators  in  Cross-Counter  Flow 


Test  names 

Test  type 

W020225B 

5 

W020228A 

6 

W020221A 

7 

W020225A 

8 

W020207B 

9 

W020530A 

10 

W020531A 

11 

W020215B 

12 

W020301A 

13 
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A. 3 Enhanced  fin  (wavy-lanced)  evaporator  in  cross-counter  flow 
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A.6  Wavy  fin  evaporator  with  non-uniform  airflow 


Table  A. 6.1 : Wavy  Fin  Evaporator  with  Non-Uniform  Airflow 


Test  names 

Test  type1 

Velocity  ratio 

W020522A 

9 

1:1 

W020523A 

9A 

1:1.5 

W020524A 

9B 

1:1.5 

W020528B 

9 

1:1 

W020528C 

9A 

1:2 

W020529A 

9B 

1:2 

])  Test  9:  uniform  airflow  with  superheat  on  all  circuits  of  5.6  °C  (10.0  °F),  9A:  same  expansion  valve  setting  as  test  9 but  with  non-uniform 

airflow  and  no  superheat  adjustment,  9B:  expansion  valves  adjusted  to  yield  5.6  °C  (10.0  °F)  superheat  on  all  circuits  with  non-uniform  airflow. 
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A. 7 Enhanced  fin  (wavy-lanced)  evaporator  with  non-uniform  airflow 


Table  A. 7.1 : Enhanced-Cut  Fin  (Wavy-Lanced)  Evaporator  with  Non-Uniform  Airflow 


Test  names 

Test  type1 

Velocity  ratio 

E020604A 

9 

1:1 

E020604B 

9A 

1:1.26 

E020605A 

9B 

1:1.26 

E020607A 

9 

1:1 

E020607B 

9A 

1:1.36 

E020610A 

9B 

1:1.36 

E02061 1 A 

9A 

1:1.62 

E020612A 

9B 

1:1.62 

E020613A 

9A 

1:1.75 

E020620A 

9B 

1:1.75 

E020621 A 

9A 

1:2.59 

E020624A 

9B 

1:2.59 

1)  Test  9:  uniform  airflow  with  superheat  on  all  circuits  of  5.6  °C  (10.0  °F),  9A:  same  expansion  valve  setting  as  test  9 but  with  non-uniform 

airflow  and  no  superheat  adjustment,  9B:  expansion  valves  adjusted  to  yield  5.6  °C  (10.0  °F)  superheat  on  all  circuits  with  non-uniform  airflow. 
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APPENDIX  B.  CAPACITY  UNCERTAINTY 


Table  B.l  lists  the  relative  uncertainty  in  the  air-side  capacity  for  two  representative  tests  at  low 
and  high  evaporator  capacity.  Two  tests  are  shown  below,  with  the  first  test  being  a typical  test 
at  a capacity  comparable  to  a majority  of  the  other  tests  for  all  coils.  The  second  test  listed  in 
Table  B.l  shows  a worst  case  test  for  COIL-W  in  parallel  flow  with  an  extremely  low  capacity. 
For  the  majority  of  tests,  the  uncertainty  in  the  evaporator  capacity  was  at  the  4 % to  5 % level. 
A complete  description  of  the  propagation  of  error  technique  used  to  calculate  uncertainty  is 
given  in  Payne  and  Domanski  (2001). 


Table  B.l:  Relative  Uncertainties  of  Two  Evaporator  1 

"ests 

Test  Name 

Coil 

Designation 

Capacity, 
kW  (Btu/h) 

Uncertainty 

Description 

Capacity, 
kW  (Btu/h) 

Percent  Uncertainty 
at  a 95  % Confidence 
Limit  on  the  Mean 

E020416B 

Enhanced- 

Cut 

7.8  (26546) 

Typical  of 
all  tests 

7.8  (26546) 

4.2 

W020311B 

Wavy 

0.90  (3078) 

Worst  case 

0.90  (3078) 

8.9 

166 


APPENDIX  C.  USER’S  INSTRUCTION  FOR  THE  EVAP-COND  VERSION 

USED  IN  THIS  STUDY 

A CD  attached  to  this  report  contains  a version  of  EVAP-COND  that  was  specifically  developed 
for  this  study.  The  following  pages  describe  how  to  install  and  use  the  model.  As  needed  for 
this  study,  this  version  of  EVAP-COND  simulates  only  evaporators  with  multiple  inlets  using 
the  option  that  solicits  refrigerant  outlet  saturation  temperatures  and  global  superheat.  This 
option  is  identified  in  the  figure  below  with  the  EVAPORATOR  OPERATING  CONDITIONS. 
The  condenser,  which  normally  is  included  in  the  EVAP-COND  package  is  not  provided  here. 

The  attached  CD  package  contains  the  following  two  files: 

EV-CD.exe  - self-extracting  file  with  all  files  needed  for 

executing  EVAP-COND. 

EVAP-COND  instructions.pdf  - file  with  visual  instructions  on  how  to  use  EVAP- 

COND.  (You  need  Version  5 of  Adobe  Reader  to 
read  this  file.).  The  instructions  are  also  included 
in  this  appendix. 

Installation  of  EVAP-COND  on  your  PC 

Execute  file  EV-CD.exe  to  expand  it  on  your  hard  drive.  You  will  be  prompted  to  select  a 
directory  where  you  want  the  program  to  reside.  When  the  installation  is  completed,  you  should 
see  EVAP-COND  directory  and  two  subdirectories  called  FLUIDS  and  MIXTURES. 
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In  the  main  directory  (EVAP-COND),  EVAP-COND.exe  is  the  interface.  EVAP5.exe  is  the 
evaporator.  Files  with  the  affix.dat  are  example  cases  of  input  data  used  in  this  study.  Files  with 
the  affix. opc  extension  contain  corresponding  operating  conditions. 


Next  step 

Refer  to  the  following  pages  or  EVAP-COND  instructions.pdf  for  further  information  about 
EVAP-COND  capabilities  and  limitations.  They  will  also  assist  you  in  your  first  evaporator 
simulation  run.  It  is  recommended  for  the  user  to  follow  the  steps  described  there  to  familiarize 
yourself  with  the  model.  Because  of  constant  upgrading  of  the  model,  the  simulation  results  you 
are  going  to  obtain  may  not  be  the  same  as  those  presented  in  EVAP-COND  instructions.pdf. 


Control  of  the  option  to  simulate  with  or  without  longitudinal  fin  conduction 

The  user  can  control  the  option  of  using  longitudinal  fin  conduction  in  a simulation  by  accessing 
file  TUNE.TXT  selecting  0 or  1 for  the  flag,  as  it  is  explained  in  the  file,  and  saving  the  file. 
TUNE.TXT  is  located  in  the  EVAP-COND  directory. 

The  following  pages  contain  general  visual  instructions  for  using  EVAP-COND  as  they  are 
presented  in  the  file  EVAP-COND  instructions.pdf  located  on  the  attached  CD.  The  option 
developed  for  this  project  is  marked  and  available  in  this  package  is  marked. 
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EVAP-  COND  INSTRUCTIONS 


EVAP-COND  is  a software  package 
that  contains  NIST’s  simulation  models 
for  a finned-tube  evaporator  (EVAP5) 
and  condenser  (COND5).  The  following 
pages  provide  basic  instructions  on  how 
to  use  this  package.  The  instructions 
include  preparation  of  input  data, 
execution  of  the  program,  and 
examination  of  simulation  results. 

Capabilities  include: 

• tube-by-tube  simulation 

• non-uniform  air  distribution 

• simulation  of  refrigerant  distribution 

• condenser  model  capable  of 
simulattions  above  the  critical  point 

•10  refrigerants  and  refrigerant  mixtures 

• REFPROP6  refrigerant  properties 


Piotr  A.  Domanski 

National  Institute  of  Standards  and  Technology 
Building  and  Fire  Research  Laboratory 
Gaithersburg,  MD,  USA 


Smart  Distributor  Version,  Aug.  5. 2002 


EVAPORATOR  REPRESENTATION 


LOADING  A FILE 


After  starting  the 
program,  click  on  the 
“open  file”  button  on 
the  power  bar  or 
select  "Open"  on  thi 
pull  down  menu 


Open  file 
EXAMPLE-5  to 
continue  this 
tutorial. 


PREPARING  A SIMULATION  RUN 
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REFRIGERANT  CIRCUIT  DESIGN 


Follow  the  steps  below  to  design  a circuitry: 

1 . place  the  pointer  on  the  tube 

2.  press  the  left  mouse  button 

3.  drag  the  pointer  to  the  next  tube 

4.  release  the  left  button 


REFRIGERANT  SELECTION 


0 P®  I £*  fiunSniatan  Satiation  Rea*:  )£**  Itfretow 


EXAMPLE-5  was  set  up  to  use 
refrigerant  inlet  pressure  and  quality  as 
input  data  (see  Operating  Conditions). 
Since  the  specified  inlet  pressure  is  for 
R22,  use  the  pull-down  menu  to  select 
R22. 
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COIL  DESIGN  DATA 


Value  1 must 

be  used  for 
this  version  of 
EVAP-COND. 


AIR  VELOCITY  PROFILE 


To  specify  velocity  at  a given  point, 
place  the  pointer  at  that  point  and 
click  with  the  left  button.  The  window 
above  shows  the  coordinates.  Up  to 
sixteen  velocity  points  may  be 
specified. 

Click  on  a given  point  to  remove  it. 

The  Velocity  Profile  window  must  be 
open  to  make  air  velocity  profile 
changes. 


to  disregard  the  previously  specified 
vol.  flow  rate  and  calculate  a new 
volumetric  flow  rate  from  the  velocity 
profile 


Recall  that  the  operating  conditions 
input  data  included  the  volumetric  flow 
rate  of  air.  Before  you  leave  this 
window  you  must  decide  between  the 
two  options: 

• to  use  the  previously  specified  vol. 
flow  rate  and  have  the  velocity  profile 
scaled  to  match  that  value 
(recommended) 


EVAPORATOR  OPERATING  CONDITIONS 

(Eight  options) 


Input  includes  estimated  refrigerant 
mass  flow  rate,  which  is  later 
iterated  by  EVAP5. 

Inlet  quality  is  imposed  at  0.20. 


Input  includes  refrigerant  mass  flow 
rate. 


This  option  includes  a refrigerant 
distributor.  Condenser  subcooling 
and  evaporator  superheat  are 
imposed  to  be  5 °C  (9  ° F). 


Input  indudes  estimated  refrigerant 
mass  flow  rate,  which  is  later 
iterated  so  the  least  superheated 
tube  has  the  specified  minimum 
superheat 

Inlet  quality  is  imposed  at  0.20. 


EXECUTION  OF  EVAP5 
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EVAP5  OPENING  WINDOW 


EVAP5  SOLICITS  REFRIGERANT 
DISTRIBUTION 


This  window  appears  only  in  this  version  of  EVA P-  COND 
for  the  simulation  option  used  for  this  study 
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SIMULATION  RESULTS 


SIMULATION  SUMMARY 


The  Simulation  Results 
Summary  window 
displays  results  in  the 
units  selected  for  the  input 
data.  Files  si.res  and 
bt^es  have  the  same 
information  in  the  SI  and 
Inch-Pound  units, 
respectively. 


SIMULATION  SUMMARY  (cont.) 

(Complete  printout  of  file  sixes) 


HOW  TO  SIMULATE  EVAPORATOR  ? 

(An  example  using  the  existing  file  EXAMPLE-5.dat) 


Run  Windows  Explorer  and  go  to  the  directory  containing  EVAP-COND.exe. 

Double-click  on  EVAP-COND.exe  to  start  the  program. 

Open  file  EXAMPLE-5.dat  to  simulate  the  evaporator.  After  the  file  is  loaded,  you  will  see  a schematic  representing  a 
side  view  of  the  evaporator.  The  red  ctrcle(s)  indicates  the  inlet  tube  to  the  evaporator  The  blue  circles  indicate  the 
outlet  tubes.  The  horizontal  line  at  the  bottom  of  the  screen  indicates  the  air  velocity  profile  at  the  evaporator  inlet. 

Review  Input  Data.  Click  on  the  Edit/Coil  Design  menu  item  to  review  the  evaporator  design  information.  You  may 
select  either  the  SI  or  British  system  of  units  for  your  input  data  and  simulation  results. 

Click  on  the  Edit/Operating  Conditions/Evaporator/iniet  pressure  and  quality  menu  item  to  review  operation 
conditions.  Note,  that  the  loaded  option  has  a mark  on  the  left-hand  side.  Since  EVAP5  simulates  performance  tube- 
by-tube  from  the  inlet  to  outlet,  the  options  that  specify  any  outlet  refrigerant  parameter  involve  iterative  calls  to  the 
option  that  specifies  refrigerant  inlet  pressure  and  quality  until  the  target  outlet  parameters  are  obtained  (e  g., 
saturation  temperature  and  superheat). 

Click  on  the  EdiEVelocity  Profile  menu  item  to  review  the  air  velocity  profile.  You  may  use  the  air  mass  flow  rate 
specified  earlier  in  the  Operating  Conditions  window  or  integrate  the  air  velocity  profile.  In  general,  the  first  option  is 
recommended  unless  very  detailed  and  accurate  local  measurements  of  the  velocity  profile  were  taken.  You  may 
change  the  air  velocity  profile  using  a mouse  by  clicking  the  left  button. 

Run  a simulation.  Click  on  the  Run  Simulation  menu  item  and  select  EVAP5  An  MS-DOS  window  will  appear  and 
will  give  you  a message  when  a simulation  run  is  successfully  completed 

Examine  local  and  global  simulation  results  EVAP5  writes  global  simulation  results  to  file  Si  res  (SI  system  of  units) 
and  BT.res  (British  system  of  units).  The  same  information  is  provided  in  the  pull-down  menu  in  the  units  selected  for 
data  input. 


HOW  TO  PREPARE  YOUR  DATA  FILE  ? 


Start  with  Edit/Coil  Design  menu  item.  Input  all  information 

Select  Edit/Operating  Conditions  menu  Item  to  Input  operating  conditions  data. 

Select  EditA/elocity  Profile  to  change  the  velocity  profile  using  a mouse  (left  button). 

Specify  refrigerant  circuitry. 

If  you  are  coding  evaporator  circuitry,  start  with  one  of  the  inlet  tubes  and  proceed  downstream.  If  you  are  coding 
condenser  circuitry,  start  with  one  of  the  outlet  tubes  and  proceed  upstream,  i.e.,  in  either  case  you  have 
to  start  from  the  side  that  is  closer  to  the  saturated  liquid  line. 

To  draw  a return  bend,  point  the  mouse  on  a tube,  press  the  left  button,  drag  the  mouse  to  the  next  tube,  and  release 

If  you  want  to  modify  a circuitry,  you  may  delete  a part  of  it  starting  from  a given  tube  end  ending  by  the  exit  tube  by 
pointing  the  mouse  on  the  given  tube  and  double-clicking  the  left  button. 

Once  a circuit  is  coded,  it  can  be  used  for  both  evaporator  and  condenser  simulations  based  on  specified  operating 
conditions. 


COMMENTS 

SUGGESTIONS 

QUESTIONS 
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